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SYNOPSIS 


This final report for NASA Grant NCC 2-711 cover reporting period June 
1992 through December 1992. The report analyzes the longitudinal and lateral 
flying qualities of Propulsive-Only Flight Control (POFC) for a Boeing 720 
aircraft model. Using Quantitative Feedback Theory (QFT), performance results 
from compensators are documented and analyzed. This report is also the first 
draft of a graduate thesis to be presented by Hwei-Lan Chou. The final thesis 
document will be presented to NASA when completed later this year. 

The latest landing metrics, related to bandwidth criteria and based on the 
Neal-Smith approach to flying qualities prediction, were used in developing the 
performance criteria for the controllers. The compensator designs were tested on 
the NASA simulator and exhibited adequate performance for piloted flight. 

There was no significant impact of QFT on the performance of POFC in either the 
longitudinal or lateral modes of flight. This was attributed to the physical limits 
of thrust available and the engine rate of response, both of which severely limited 
the available bandwidth of the closed-loop system. 


viii 


Abstract 


Through throttle manipulations, engine thrust can be used for emergency 
flight control for multi-engine aircraft. Previous study by NASA Dryden has 
shown the use of throttles for emergency flight control to be very difficult. In 
general, manual fly-by-throttle is extremely difficult - with landing almost 
impossible, but control augmentation makes runway landings feasible. Flight 
path control using throttles-only to achieve safe emergency landing for a large jet 
transport airplane, Boeing 720, was investigated using Quantitative Feedback 
Theory (QFT). Results were compared to an augmented control developed in a 
previous simulation study. The control augmentation corrected the 
unsatisfactory open-loop characteristics by increasing system bandwidth and 
damping, but increasing the control bandwidth substantially proved very 
difficult. The augmented pitch control is robust under no or moderate 
turbulence. The augmented roll control is sensitive to configuration changes. 
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1. Introduction 


Through throttle manipulations, engine thrust was found useful in providing 
some controllability for multi-engine aircraft in emergency situations with severe or 
complete flight control system failures (such as hydraulic system failures). Aircraft 
flight control systems are extremely reliable. Current generation aircraft utilize multiple 
and independent control surfaces, hydraulics, sensors, control computers, and control 
cables to achieve a high level of control system redundancy and reliability. Although 
rare, severe flight control system failures do occur. 

NASA Dryden has studied the use of throttles for emergency flight control for a 
range of airplanes 1 - 5 . Many multi-engine airplanes exhibited some degree of useful 
control capability with the throttles. In general, flying an aircraft in manual mode using 
throttles-only requires a tremendous pilot workload and landing is considered 
extremely difficult to almost impossible. Control augmentation, using feedbacks and 
direct coupling of the throttle command to stick/ thumbwheel motion, has greatly 
improved flying qualities, and ground simulation landings can be achieved. 

The primary aim of this current study on Throttles-Only Flight Control (TOFC) is 
to develop an augmented flight path control using throttles-only to achieve safe 
emergency landings. Application of TOFC on a large four-engine jet transport airplane, 
Boeing 720 (B-720) (Figure 1), is investigated. An augmented B-720 TOFC, developed 
and implemented on a high fidelity B-720 flight simulator (Figure 2) by NASA Dryden 2 , 
had obtained good pilot rating by increasing the control bandwidth and the phugoid 
and Dutch-roll damping 2 . 

This report presents an alternative control design technique based on 
Quantitative Feedback Theory (QFT) to further improve the Dutch-roll damping and to 
increase the control bandwidth for better handling qualities. The control design uses a 
linearized B-720 model derived from perturbations of the full non-linear equations of 
motion about trim at an approach and landing flight condition. 

A robust controller is highly desirable for systems with plant parameter 
uncertainty (such as an aircraft undergoing configuration changes). The QFT 
technique 6 - 9 was chosen because it allowed designers to specify a desired close-loop 
response and a performance specification, and then built a controller to meet the 
specification. Most of all, the technique can incorporate plant parameter uncertainty 
and plant disturbances into the control system design by converting them into design 
constraints and then design a controller to have the system satisfy the imposed 
constraints. The controller thus designed guarantees robust performance over full 
range of the plant uncertainty while keeping the disturbance effect to the system 
minimum. 

The desired performance specification may not always be achieved within the 
given control actuation and rate limits. However, the transparency of the QFT 
technique throughout the design process preserves many of the insights which are lost 
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in several of the modem control techniques and thus provides control designers with 
valuable information about the system under investigation. QFT also provides a 
quantitative relationship between the amount of uncertainty and feedback (i.e. the 
magnitude of feedback is determined in proportion to the amount of uncertainty, 
therefore, reduces the possibility of overdesign). 

In this report, the strategy of flight control using throttles-only is introduced. 
The fidelity of the linear B-720 model is examined. An overview of QFT with step-by- 
step procedures is provided, and its application on the design of an augmented flight 
path control using throttles-only for approach and landing of B-720 is presented in a 
summary fashion. Control design results using QFT are compared to the augmented 
control developed in a previous simulation study. 


Nomenclature, Abbreviations, and Acronyms 

(a) short form of (s+ a) 

c.g. center of gravity 

C, p non dimensional yaw-roll coupling derivative 

C mu non dimensional velocity-pitch coupling derivative 

D. R. Dutch-roll 

G"" 1 transfer functions 

K q pitch rate feedback gain 

K r flight path angle feedback gain 

Kp sideslip angle feedback gain 

K # bank angle feedback gain 

q pitch rate (deg / sec) 

QFT Quantitative Feedback Theory 

s.p. short period 

TOFC Throttle-Only Flight Control 

z thrust (lbs) 

6 engine rpm 

6tc stick input (full deflection=l unit) 

y flight path angle (deg) 

0 pitch angle (deg) 

p angle of sideslip (deg) 

<p bank angle (deg) 

to i natural frequency 

£ damping ratio 

[£,w n ] short form for S 2 + 2£a> n S + a ) a 2 
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2. Strategy of Throttles-Only Flight Control 


The propulsion system of a multi-engine aircraft can be used for heading and 
flight path control. Differential throttles is applied to control roll through yaw, and 
symmetric throttles is applied to control pitch. Speed control by throttles becomes 
ineffective when control systems fail. Other means may be used to change the airplane 
speed as described in below. Throttles are coupled to stick/thumbwheel for easier and 
more conventional control handling. 

2.1 Pitch Control 

Symmetric throttles induces a phugoid mode and a speed change, which in turn 
generates a pitching moment change through speed stability effects-C mu . This is the 
primary source of pitch control. Pitch control may also be generated by other factors 
such as pitching moment change due to thrust line offset, flight path angle change due 
to the vertical component of thrust, and an instant pitching moment change generated 
by engines mounted at different vertical levels, as in the case of B-720. 

2.2 Yaw-Roll Control 

Differential thrust generates sideslip, which in turn generates rolling moment 
changes through wing dihedral and sweep effect-C^. Roll is controlled by applying 
differential throttles to achieve the desired bank angle and thus to make turns and 
heading changes. 

2.3 Speed Control 

Retrimming speed by the use of throttles becomes ineffective when primary 
control surfaces are locked due to control systems failure. When control system failure 
occurs at speeds other than landing speed, retrimming to an acceptable landing speed 
maybe accomplished by using other techniques such as lowering flaps (assuming the 
electrically controlled flaps are operative), extending landing gears, moving eg. aft, 
varying stabilizer deflection, or varying the speed between the low and high mounted 
engines. 

2.4 Couple Throttle Command to Stick/Thumbwheel Motion 

Direct coupling of the throttle command to stick/thumbwheel motion has eased 
the pilot's handling of control. The airplane can be controlled in a conventional fashion, 
such as pitch up with stick forward or pitch down with stick aft. 
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3. B-720 Linear Model 


The B-720 linear model is derived from perturbations of the full nonlinear equation 
of motion about trim and is completely decoupled in longitudinal and lateral dynamics. 
All control states and inputs are perturbed independently at the steady state of a 
desired trim condition. The inputs are thrust from each engine. There are four 
configurations given for the study of approach and landing of B-720 TOFC. The state- 
space representation of the linear model of these configurations are listed in Appendix 
A. Of the four configurations shown, configuration 1 is the nominal configuration for 
baseline design. 

The fidelity of the linear model was examined by comparing the open-loop 
response of the linear model with the nonlinear model as shown in Figures 3 and 4. The 
longitudinal response of the linear model was about 30% less in magnitude than the 
nonlinear model, and was, therefore, modified by a correction factor of 1.3 (Figure 3). 
The longitudinal linear model after modification closely portraits the longitudinal 
nonlinear model (Figure 5). Linear design analysis utilized a computer control package 
"Program CC" to assist the design 10 . 

The lateral response depicted in Figure 4 shows that the linear model would closely 
follow the nonlinear model as long as the small perturbation assumption is not violated, 
i.e. the command input should be of a small magnitude and a short duration. Figures 3 
and 4 also illustrate that for the nonlinear model, a flight path angle command would 
induce little coupling in roll /yaw, while a bank angle command would induce 
pronounced pitch coupling. Coupling between longitudinal and lateral modes is 
completely absent for the linear model. 


4. Engines and Bare Airframe System Analysis 
4.1 Engines 

Spool-up and spool-down engine dynamics for the B-720 engine are shown in 
Figure 6a 12 . The empirical transfer function developed is given in short form notation 
by 


QZ(lbs) _ 
a Tt(*) 


275 

(•55)(5) 


This equation is illustrated in Figure 6b over low frequency ranges up to 1.0 
rad /sec. 
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4.2 Bare Airframe 


It is apparent from the engine Bode diagram in Figure 6b illustrates that severe 
bandwidth attenuation would occur beyond frequencies of 1 rad /sec. Therefore, it may 
not be possible to increase the closed -loop bandwidth beyond 1 rad / sec within the 
range of available thrust. This can be seen in the pitch rate " o[' to thrust "z" transfer 

function, G^ (de8/,ec) (refer to Appendix A), of the bare airframe shown in Figure 7. The 
full-order transfer function Gf deg/ ‘ ec) shows that 80 dB of gain must be added to yield a 

crossover frequency beyond 1 rad / sec. This corresponds to 10,000 lbs of full thrust 
from each engine, which is not practical for approach and landing. 

A low order fit to Gf deg/ “ c) is also depicted in Figure 7 and is very accurate near 

the phugoid frequency. Piloted flight of the unaugmented aircraft was consistently a 
level 3 2 . The main difficulties were the lightly damped phugoid and the low 
bandwidth throttle control. 

The accuracy of the low order fit near the phugoid frequency means that, to a 
first order approximation 12 , the phugoid frequency and damping are found from the 
following equation: 


2£c0 o - -X u + 


M u (X a -g) 


M 

-g(Z u -^Z a ) 


w„ 2 - 




U„ 


Conventional transport aircraft can be shown to be roughly proportional to M u . 

It should be strongly noted here for the classic case of M u =0 and for negative 
values of M u (Mach tuck) that the aircraft cannot be practically flown with throttle alone 
unless rotational control in pitch is added and difficulties will also be encountered as 
M a becomes small (aft c.g. location). Both of these cases require the addition of an 
effective rotational controller about the pitch axis. This may be achieved by using 
differential inboard and outboard thrust, provided the inboard engines are a different 
distance from the aircraft xy-plane than the outboard engines. These configuration 
characteristics determine the innate capability for throttles-only piloted control. 
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5. Overview of Quantitative Feedback Theoiy 

QFT is a frequency domain control technique that uses a fairly straightforward 
and transparent design approach 6 - 9 . To apply QFT, systems are usually modeled in a 
unit feedback form (Figure 8) where all blocks may present scalar (SISO) or matrix 
(MIMO) system transfer functions. For MIMO systems, a m x m MEMO system can be 
converted into am 2 - equivalent multi-input single-output (MISO) loops (Figure 9). 

QFT techniques allow designers to specify a desired performance specification with 
performance tolerance and then incorporate the tolerance with the plant uncertainty 
and system disturbances to form the design constraints: the performance bounds and 
the U contour (Figures 10, 11 and 12). 

The design constraints are then placed on a Nichols Chart together with the 
nominal plant transfer function, P G . A controller will be selected to reshape P Q to form 
Lq (the nominal open-loop transfer function) as to have Lq satisfying all the design 
constraints of performance bounds and U contour (Figure 13). By having L 0 satisfy all 
the design constraints, if possible within the given control actuation and rate limits, the 
system is guaranteed robust over the full range of plant uncertainty. However, the 
system may not completely meet the performance specification (Figure 14). A prefilter 
is usually required to further reshape the system to fully meet the specification. The 
prefilter design is implemented on a Bode plot. 

The basic design procedures of the QFT technique for minimum phase systems 
are accomplished by the following four steps: 

1) Model the system in a unit feedback form to apply QFT. A m x m MIMO system can 
be converted into a m 2 -equivalent MISO system and the coupling between loops can 
be considered disturbance input (Figures 8 and 9). 

2) Specify the desired close-loop frequency response performance specification. Figure 
10 shows the construction of a desired close-loop performance specification with an 
upper bound, Bu ; a lower bound, Bl; a tolerance, and a maximum peak 
magnitude, M m . 

The tolerance, <5 H , is specified to obtain robust performance, and the maximum 
M m is specified to obtain a desired system damping. The upper bound is generally 
synthesized by an underdamped second order close-loop transfer function (T.F.), 
T u (s) and the lower bound by an overdamped close-loop T.F., T L (s) with figures of 
merit such as settling time, rise time, peak overshoot or damping ratio, and natural 
frequency, etc. 

A desired disturbance performance specification (Figure 11) needs only an upper 
bound to confine the disturbances. The objective of the technique is to design a 
controller such that the variation of the response due to plant uncertainty lies within 
the specified boundaries and the effect of disturbance is minimized, that is to have: 
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Bl(cd) £ |TrQ o>)| s Bu(w) 

|TdiOco^|MdO(o)| 

|T D2(j to)j £ |Mo(j (o)| 

3) Convert the performance tolerance, <5 n , and the maximum onto Nichols Chart to 
form the design constraints: the performance bounds, Bo(jG>i), and the U contour. 

i) Performance bounds are curves on the Nichols Chart that are determined by 
matching the magnitude of the range of plant uncertainty with the magnitude of 
the performance tolerance, <5„. Therefore, satisfying this constraint guarantees 
the variation of the system response due to plant uncertainties will be no greater 
than <5 r . There is a performance bound for each frequency. 

ii) On the Nichols Chart, the U contour is a M-circle that has the magnitude of M m , 
with part of the circle stretched for uncertainty at high frequencies (same as the 
length V shown in Figure 12). By having the open-loop response not penetrating 
the U contour, the system's damping will be guaranteed no less than the 
damping correlating to Mm- The construction of a U contour is shown in Figure 
12 . 

4) Reshape the nominal plant transfer function, Po. Gain /pole/ zero compensation 
may be placed on Po to reshape it to satisfy the design constraints. After reshaping, 
Po becomes Lo, and the compensation chosen forms the controller, G 0 as can be 
depicted from the relationship: 

L 0 = P 0 * G c 

To satisfy the design constraints, Lq should not penetrate the U contour, while each 

frequency co, on L 0 should be kept on and above its corresponding Boflcoj). The U 
contour , the performance bounds and the optimal Lo of an example problem are 
shown in Figure 13. 

5) After reshaping, the system is guaranteed robust over the full range of plant 

uncertainty, i.e. &r(jwi) s ^R(j w 0 (Figure 14). However, the system may not have 
met the performance specifications completely. A prefilter is usually required to 
further reshape the system to fully meet the specification. 
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6. Quantitative Feedback Theory Control Design 

A QFT computer control package was used to assist the QFT design 8 . The 
program is to be used for minimum-phase plants only, i.e., the plants should have no 
zeros in the right half s-plane, therefore, only the gain curve of the desired performance 
will be specified and satisfied. For nonminimum-phase plant, the phase of the desired 
close-loop performance shall also be specified and satisfied. 

6.1 System Modeling 

The block diagrams of flight-path-angle control and bank-angle control are 
presented in Figures 15 and 16. The inner pitch-rate (q) loop and sideslip-angle ((5) loop 
were first closed with G*| c =1, K q =60 and G^ c =10, K^=4, respectively, which were the 

heuristic settings chosen by investigating the properties of the inner loop. Tables 1 and 
2 summarize the investigation: 


Table 1. Investigation of Longitudinal Feedback Parameter 


Feedback 

Parameter 

Phugoid Mode j 


CD n 

q 

Increase 

(Require high gain) 

No change 

Y 

No change 

Increase 


Table 2. Investigation of Lateral Feedback Parameter 


Feedback 

Parameter 

Lateral Phugoid 

Dutch Roll 

% 

0) n 



P 

Increase 

Increase 

Decrease 


r 

Small 

increase 

No change 

Small 

increase 

No change 

0 

Small 

decrease 

Small 

increase 

Increase 

No change 

<p 

Increase 

Increase 

Small 

increase 

No change 


To apply QFT, with the inner loop closed, the outer loops are rearranged in a unit 
feedback form as shown in Figures 17 and 18. 

































6.2 Performance Specification 

To obtain good handling qualities, the close-loop response for each of the y- and 
<p- loops, which are also the pilot control open-loops, should have the following 
characteristics: 

1) A bandwidth, (D^b, - 2 rad /sec for landing of a transport aircraft 11 . 

2) A k/s gain curve slope (-20 dB /decade) around the crossover frequency, w c u . 

A desired close-loop specification was synthesized based upon these two 
requirements, and is shown in Figure 19, which has a k/ s slope near C0c — 1*5 rad /sec 
(with (Oojb - 2 rad / sec) and a comer frequency, co CO mer - 0-8 rad/ sec. 

The desired close-loop specification is synthesized in the following four steps: 

1) Synthesize the initial Bu and B L . Bu is usually modeled by an underdamped 
simple second order close-loop T.F., 

T u (s)= 2 - 7 , 

S + 2gwr„S + tu a 

while B l modeled by an overdamped simple second order close-loop T.F., 


T, (S)= , wherecr, and a 2 s tj - 

(8 + ^XS+aa) 

With a desired performance specification of g=.6 and cr n =.8 rad /sec, this 
yields: 


.385 


T u( s )“ s 2 + 96s+ 64 311(1 Tl ^“ (S+.55XS+.7) 

2) Add a pole to T L (S) to widen the <5„ at high frequencies. This yields: 


T l (S)- 


1 77 

(S+.55XS+.7XS + 2) 


This is required by the Bode derived theorem which states that |l m Sjdc7 - 0 , 

i.e. the reduction in sensitivity Sj at the lower frequencies must be 
compensated by an increase in sensitivity at the higher frequencies. 
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3) Add a zero at 1 rad /sec to increase the gain slope from -40 dB/decade to -20 
dB /decade. This yields: 


T„(s) 


: 64(S + 1) and! 


L (S)' 


.77(S + 1) 


S 2 +.96S+.64 (S+.55XS+.7XS + 2) 

4) Raise the whole synthesized gain curve until a cOo*, « 2 rad /sec is obtained. 


The magnitude of B u , B!,and <5 n at each frequency can thus be obtained and are 
shown in Table 3. 


Table 3. QFT performance specification 


Frequency(rad /sec) 

0.1 

0.3 

0.5 

mm 

1.0 

2.0 

5.0 

Bu(dB) 

17.0 

17.0 

17.3 

16.0 

13.0 

2.0 

-13.0 

BL(dB) 

16.8 

15.0 

12.3 

9.9 

4.6 

-7.1 

-23.0 

<5„ (dB) 

0.2 

2 

5 

Em 

Em 

9.1 

15.0 


6.3 Airplane Parameter Uncertainty 

Four configurations are provided for the study of approach and landing of B-720 
throttles-only flight control. The flight condition of these configurations are 
summarized in Table 4. Configuration 1 is the nominal configuration for baseline 
design. 


Table 4. Flight Configurations for B-720 A 

pproach and Landing 

Config. 

Number 

Weight 

(lbs) 

Altitude 
(Ft MSL) 

Airspeed 

(Knots) 

Flaps 

(%) 

Gear 

up /down 

1 

140,000 

4,000 

160 

0 



2 

140,000 

4,000 

145 

30 



3 

160,000 

4,000 

175 

0 

ilE 

4 

140,000 

4,000 

155 

30 

_HE 
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A plant transfer function with parameter uncertainty is usually described in a 
maximum and minimum format in order to form the plant uncertainty template, which 
will then be used to determine the performance bounds constraint. An example of a 
plant wuth parameter uncertainty described in a maximum and minimum format is 
shown below: 


Example 


For a plant transfer function 


G(S)« 


Ka 

(s + a)’ 


where the parameter variations are: 1 < k < 10 and 1 < a < 10 


then, 

G(S) min. - — — and 

S + l 


G(S) max. - 


100 

S + 100 


For y- and tp- feedback loops, the minimum and maximum values of the transfer 
functions G^ and G£ , determined from the four given configurations, are shown 
below: 


For y-feedback loop: 


The G^ ( ^ g) of the nominal configuration(config. 1) is: 


G£3L)«*nfig. 1 = 


.0l(.203)[.37, 3.01 


(.562) [.624, .111] [.441, L57] (5.25) 


and the min. and max. G^^ are: 

cr(*g) . .0053 (.162) [.35,3.01] 

(.40) [.42, L48] [.66.. 01] (5.19) 

-01 (.28)[.46, 3.43] 

‘ (,58)[.45, 157][.92, .14] (5.24) 
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For 0 - feedback loop: 


The G* ( *®? g) of the nominal configuration(config. 1) is: 

nominal = •. 09 [- 4 7 .^g] 

V<*.) nommal (. 98 ) [, 81 , . 15 ] [.26 , L07] (5.02) 

and the min. and max. of G*^®^ are: 


min = 


G Ji d S,) max -= 


. 06 [. 45 , 3.65] 

(.98) [.60, .15] [.24 , .93] (5.0l) 

.09 [.61 , 4.33] 

(1.03) [1.0 , .20] [.29 , L09] (5.02) 


The QFT control package, used to assist the design, allows the designer to input 
plant parameter variations by entering the transfer function's maximum and minimum 
values for gain, first order poles and zeros, and second order poles and zeros. The 
program forms the plant uncertainty template with the given maximum and minimum 
values, then uses its CAD capability to graphically determine the performance bounds 
required for the design. There are tradeoffs between plant parameter uncertainty and 
system performances. The wider the spread of the parameter uncertainty, the more 
restricted the constraints; consequently more compensation is required. Therefore, the 
performance specification may need to be relaxed when there is not enough control 
power to provide all the compensation that is required. 


6.4 Controller and Prefilter Design 

Pole/zero/gain compensation maybe required to reshape the plant transfer 
functions of the y- and tp- feedback loops, G^ and G^, and to satisfy performance 

bounds and U contour constraints. On a Nichols Chart, adding a gain will raise the 
transfer function curve, while a zero will bend the curve to the right, and a pole will 
bend the curve to the left. The compensation chosen forms the controller, G c . After 
reshaping G^ and G^ become, respectively, L r er and (the open-loop transfer 

functions of the y- and <p-loops), where L r er - G®‘ r " * and - G* 1 " * G^ . Each 
frequency coj on and should be placed on and above its corresponding 
performance bounds, Bo(ja)j), to assure robust performance. In addition, L^ r and 
must not penetrate the U contour in order to obtain the desired damping. 

Longitudinal Flight-Path- Angle Control: Transfer function Gj^ and its performance 

bounds, Bo(ja)i), and U contour are displayed on a Nichols Chart in Figure 21. All 
frequency points on G^ n are below their corresponding Bo(ja>i), hence reshaping is 
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required (Figure 21). A pure gain compensator, G C =G® 1 " =16, raises the curve to just 
touching the U contour (Figure 22). Several lead compensators were tried to further 
reshape the G£ ( to satisfy all BoQo),) while not penetrating the U contour. The lead 
compensators tried had increased the bandwidth and robustness; however, they also 

reduces the output y to a very small value(e.g., an output y=2 degree for a full stick 
input). Therefore, only the pure gain of 16 is chosen as the compensator, this left the 

Bo(j<Oi) unsatisfied. The frequency response of the close-loop transfer function, 

T r r ta (where = L r er / (1 + L y er ) - (Gj; *G^)/(1 + Gj; * G^ ), is shown in Figure 23. It can 
be seen in Figure 23 that d T (the spread between T ma x and T m in) had exceeded the <5 R 
over the frequency range 0.1 to 0.7 rad / sec as a result of L^ y not satisfying the 
performance bounds over that frequency range. To have any frequency, <o if on L' r 
higher than its corresponding Bo(ju>i) will result in <5 T (ja>,) z <5 B (j<w,), while lower than 

Bo(jwi) will result in <5 T (j(u ( ) se <5 B (ja>|). As can be depicted in figure 23, further 
modification is required to fully meet the prescribed specification. A prefilter of a pure 
gain of 6.5 proved most effective in increasing the bandwidth and met the prescribed 
specification. The frequency response after the prefilter is applied is shown in Figure 
24. 

Lateral Bank-Angle Control: Transfer function G^ and its performance bounds, 

Bo(jtOi), and U contour are displayed on a Nichols Chart in Figure 25. Notice that G^ is 
not only below all performance bounds Bo(jtOi) but it also penetrates the U contour. 
Therefore, more than just a pure gain is required to reshape G^ . A controller, 

G^J =(s+.15)/(s+1.5), was added to G^ to reshape it and prevent it from penetrating the 

U contour, but it was not successful in satisfying all of the Bo(joDj). After reshaping, 
is shown on a Nichols Chart in Figure 26. The frequency plot of the close-loop transfer 
function, T£ , where = L*, / (1 + L^) - (Gj* * GJ^ )/(Gjj * Gj^ ), with no prefilter 
applied yet is shown in Figure 27. A lead compensator of (S+l)/(S+2) is added to 
haunch up the severely deteriorated curve at frequency over 1 rad /sec and to increase 
the phase margin. A lag compensator of (S+0.25)/(S+0.15) is added to steepen the gain 
curve at low frequencies and to provide a smoother k/ s curve for good pilot handling 
qualities. The close-loop response after adding the prefilter is shown in Figure 28 and 
the prefilter selected is 15(S+0.25)(S+l)/((S+0.15)(S+2)). 
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7. Results and Discussion 

The objective of this study is to improve the handling qualities for the approach 
and landing of B-720 TOFC by increasing the control bandwidth and the light Dutch- 
roll damping. The control bandwidth of TOFC depends primarily on the engine 
response to throttle command, and on the propulsion-induced low-frequency speed and 
dihedral stability effects, which are configuration-dependent, thus are fixed and 
unalterable. Therefore, the control engineer's only tools are compensation and 
feedback. 

To improve control bandwidth is very difficult as can be depicted from the pitch 
rate to thrust bode in Figure 7. Full thrust of 10,000 lbs from each engine is required to 
yield a crossover frequency just beyond 1 rad /sec. This clearly shows how control 
bandwidth is limited by the control power (the engine) available. 

For flight-path-angle control, pitch rate feedback was effective in increasing 
phugoid damping while y feedback was effective in improving frequency of the 
phugoid mode. For bank-angle control, {3 feedback was found most effective in 

increasing Dutch-roll damping while cp feedback is crucial to lateral phugoid damping. 
Yaw rate feedback, which is effective in damping Dutch-roll when rudder power is 
available, helps Dutch-roll damping and the lateral phugoid damping very little. Tables 
5 and 6 compare the dynamic modes of the bare airframe with the dynamic modes of 
previous simulation designs, of QFT design and of heuristic design(heuristic 

compensation will be discussed later on page 19.) Transfer functions of y to stick and cp 
to stick for all the four configurations are listed in Appendix A. 


Table 5. Longitudinal Mode Comparison 



Density 

Phugoid 

Short 

Period 

Engine 

Pre- 

filter 

G S| " 

7 


n 

a 

Bare 

Airframe 

(1.4E-6) 

(.04,. 13) 

(.65,1.4) 

(.55X5) 

- 

■ 

■ 

■ 

■ 

Simulation 

Augmented 

Control 

(4.7E-6) 

(.52,.24) 

(.52,1.5) 

(-4X5.2) 

10 

B 

10 

B 

4 

QFT 

Augmented 

Control 

(3.4E-6) 

(.62,. 32) 

(.46,1.6) 

(.3)(5.2) 

6.5 

16 

fl 

B 

60 
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Table 6. Lateral Mode Comparison 




Dutch 

Roll 

RoU 

Engine 

Pre-filter 

Gf" 


H 

H 

K, 

B 

Bare 

Airframe 

(1.1E-4) 

(.12,. 99) 

(1) 

(.55)(5) 

- 

■ 

■ 

1 

B 

B 

fl 

Simulation 

Augmented 

Control 

(.73,. 35) 

(.15,. 99) 

(1) 

(5) 

40 

1 

i 

.5 

.5 

B 

B 

QFT 

Augmented 

Control 

(.39) 

(.29,1.0) 

(1.5) 

(.45)(5) 


(.15) 

■ 

1 

1 

4 

1 

(.15X2) 

(1.5) 

Heuristic 

Augmented 

Control 

(.75,.28) 

(.22,1.0) 

(-9) 

(5) 

40 

1 


.15 

.5 

3 

B 


For longitudinal control, pure gain compensation was used. Since the short 
period mode has a frequency around 1.5 rad /sec (which was beyond the frequency that 
throttles can control) the primary concern was to increase phugoid damping and 
frequency. The phugoid damping and frequency increased from 0.52 to 0.62 and from 
0.24 rad / sec to 0.32 rad / sec, respectively. This increase of response frequency can also 
be depicted from the flight path angle response shown in Figure 29. 

For lateral control, pole /zero compensation was used. The Dutch-roll damping 
was almost doubled, from 0.15 to 0.29. The simulation augmented control has a lateral 
phugoid mode [0.73, 0.35] which combines the spiral and the slow engine mode. This 
was replaced using the QFT design with two real root modes, (0.39) and (0.45), both 
with higher frequencies, therefore faster responses. The comparison of the responses is 
shown in Figure 30. 

All plots in Figures 21 through 32 were obtained from nonlinear simulation runs 
at approach and landing conditions with major control surfaces (ailerons, elevator sand 
rudders) locked while the electrical and mechanical systems, and the landing gear 
remained operative. Figure 33 through 35 shows the throttle response, and flight-path- 
angle and bank-angle tracking response to full stick deflection. 

Turbulence Response: The response of the flight-path-angle control under 
intermediate turbulence is presented in Figure 31. Because of gust randomness, more 
than one simulation run was made to examine the tracking integrity under turbulence. 

The bank angle tracking by QFT design does not perform well under 
intermediate turbulence (Figure 32). This could be caused by the larger Kp gain (Kp=4) 
being used in the (i -feedback loop by QFT design, while a Kp=l is used for the 
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simulation designed control. This larger Kp multiplies the gust-induced-sideslip four 
times before it was fed back to the airplane. This had a dramatic effect on the bank 
angle output due to 0.2 degree of sideslip angle would generate approximately 10 

degrees of bank angle, owing to the large Qp of B-720. The sideslip angle (P ) feedback 
is the only parameter that can effectively increase Dutch-roll damping for B-720 TOFC 
A compromise seems necessary between lateral bank angle tracking and Dutch-roll 
damping. 

Good Dutch-roll damping is associated with disturbance excitation of the lateral 
phugoid mode and results in poor tracking of bank angle. Lowering the feedback gain 
reduces Dutch-roll damping but also decreases the sensitivity of bank angle to 
disturbances, and thus makes the lateral phugoid mode less troublesome to the pilot 
when flying in turbulence. 

During the investigation, it was found that the p being fed back into the B-720 
simulator was the p at the c.g. instead of the p at the nose boom. The nose boom p is 

actually measured and fed back into a real airplane. The nose boom p was then 
modeled into the B-720 simulator and the results of the bank angle tracking under 
turbulence were fairly good. Figure 36 shows the bank-angle tracking under turbulence 

due to c.g. p, while Figure 37 shows two runs of bank-angle tracking under turbulence 

due to nose boom p. The p at the nose boom has two more terms caused by lateral and 
longitudinal offsets from the c.g. The dominant term is a function of roll rate. When 
this extra term was active in the feedback loop, lateral performance improved. 

However, there was some question concerning the correctness of the sign of the yaw 
rate term as implemented in the simulation. The effectiveness of yaw rate feedback for 
TOFC need to be further investigated. 

A augmented control scheme heuristically determined that feedback p,r, p and 

<p is investigated. The compensation of this control scheme is shown in Table 6. The 
yaw rate feedback is included in the control to improve bankangle tracking, the roll rate 
and bank angle feedbacks are included to increase the damping and frequency of the 
lateral phugoid mode while the p feedback to increase the Dutch-roll damping. The 
heuristic augmented control has a slightly slower response speed ( to n =0.28 < con=0.35 of 
simulation augmented control), but a higher Dutch-roll damping (£=.22 > £=.147 of 
simulation augmented control) which has successfully damped the Dutch-roll oscillation 
as can be seen in Figure 35. 

The system response to configuration variations for y-control and for <p-control 
are shown in Figures 38 and 39, respectively. The robustness of the flight-path-angle 
control is improved by QFT as shown in Figure 38. The Dutch-roll oscillation in the 
original simulation compensation is taken out by QFT compensation; however, the 
tracking was not improved (Figure 39). Among the three augmented control developed, 
the heuristically determined augmented control presented the best robust performance 
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(with fairly good bank-angle tracking and no Dutch-roll oscillation.) Due to time 
constraints, the heuristic augmented control was not tested on the real-time nonlinear 
B-720 simulator. 


8. Conclusions and Recommendations 

Studies by NASA Dryden has shown that throttles can be used for emergency 
flight control. Manual fly-by-throttle is extremely difficult with landing almost 
impossible, but with control augmentation, runway landing is feasible. 

Flight path control design using throttles-only to achieve safe emergency landing 
for a transport airplane, Boeing 720, was investigated. Augmented throttles-only flight 
path control built in a previous simulation study has made successful simulation 
landings. However, it showed light Dutch-roll damping and low control bandwidth. 

To increase the control bandwidth substantially proved very difficult. Differential 
throttles to engines mounted at different vertical levels to generate an instant pitching 
moment may be an effective way to increase the control bandwidth. 

For throttles-only pitch control using QFT, the control bandwidth , tracking and 
control robustness were improved by QFT. For bank angle control, QFT has improved 
the Dutch-roll oscillation. However, the lateral phugoid becomes sensitive to 
configuration changes. A compromise is required between Dutch-roll and lateral 
phugoid damping given limited control power. Further investigation of the effects of 
yaw rate feedback is recommended. 
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Figure 1. Boeing-720 
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Figure 2. Boeing-720 simulation cockpit 
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Figure 3. Comparison of the longitudinal open-loop response of the B-720 

linear and nonlinear model, 20% step throttle cmd, nominal configuration 
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Figure 4. 
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Comparison of the lateral open-loop response of the B-720 linear model and 
nonlinear model, 5% differential throttle, nominal configuration 
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Figure 5. Comparison of the longitudinal open-loop response of the modified B-720 linear 
and nonlinear model, nominal configuration 
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Figure 6. Engine spooling block diagrams and Engine <5 to <5 Tc 
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Figure 7. Pitch rate to thrust bode 
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Di Da 

where 

R ; Command Step Input 
D : Disturbance Inputs 
C : output 
P : Rant 
Qc : Compensator 
F ; Prefilter 

The open loop transmission function, L, Is defined as: 

L=Qc*P and LmL = Lm<3c + LmP 

The close loop transfer functions are: 


Tracking: 

Tr - 

FL 

and Lm Tr = Lm F + Lm — 

(Di=D2=0) 


1 + L 

1+ L 

Disturbance: 

Tdi = 

P 


(R=D2=0) 


1 + L 


Disturbance: 

TD2 * 

1 


(R=D1=0) 


1 + L 



Figure 8. QFT unit feedback control structure 






Figure 9. Equivalent MISO systems of a 3x3 MIMO system 


27 







Figure 10. Typical close loop tracking specification 



Figure 1 1 . Typical disturbance rejection specification 
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Figure 14. A prefilter is required to fullfill the performance specification 



Figure 1 5. Flight-path-angle control block diagram 
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Figure 16. Bank-angle control block diagram 
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Figure 17. Flight-path-angle control with inner q-loop closed 
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Bank Angle Loop Sideslip Angle 



Figure 18. Bank-angle control with inner 0-loop closed 



Figure 19. Performance specification in time domain for approach and landing of 

B-720TOFC 
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Figure 20. 



Performance specification in frequency domain for approach and landing of B- 
720 TOFC 
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Figure 21 . Transfer function GJ , its performance bounds B(jw), and U contour on Nichols 
Chart 
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Figure 22. Open-loop transfer function , on Nichols Chart 
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Figure 23. Frequency plot of the close-loop transfer function T* with no prefilter 
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Figure 26. Open-loop transfer function, 


on Nichols Chart 
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Figure 28. Frequency plot of the close-loop transfer function with prefilter 
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Figure 29. B-720 augmented control, step flight-path-angle response with no turbulence, 

nominal configuration 


42 


- - : Simulation compensation 
— : QFT compensation 


Lateral Response 

5 degree bank angle cmd 





Time (sec) 


Angle of 
sideslip 
(deg) 


0.0 -- 


•vo 


o 


Pitch rate 
(deg/sec) 00 _ 



i 

•25 *r 


D 20 *0 .50 

Tlme(sec) 



Figure 30. B-720 augmentedcontrol, step bank-angle response with no turbulence, nominal 

configuration 
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b) QFT Compensation 

Figure 31. B-720 augmented control, step flight-path-angle response with turbulance, 

nominal configuration 
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Figure 32. B-720 augmented control, step bank-angle response with turbulance, nominal 

configuration 
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Figure 33. Right-path-angle and throttle response to full-forward stick deflection, nominal 
configuration 



Figure 34. Bank-angle and throttle response to full -right stick deflection - QFT 
compensation, nominal configuration 
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Figure 35. Response to full-right stick deflection - Heuristic compensation, nominal 
configuration 
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re 36. Bank-angle response under turbulance due to c.g. p 
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Figure 37. Bank-angle response under turbulence due to nose boom p 
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Figure 38. Robustness of the flight-path-angle control 
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b) QFT Compensation 
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Figure 39. Robustness of the bank-angle control 
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APPENDIX A: B-720 CONFIGURATIONS 


The B-720 piloted simulation can be represented by the following block 
diagram: 




The "AIRCRAFT" in the box above represents both the engine and the 
bare airframe dynamics. The engine is approximated by a transfer function, 
G en& and the bare airframe dynamics are represented mathematically by a 
single quadruple, Pa/c, shown as follows: 



Longitudinal Dynamics 


Throttle command 
% rpm 
<5 t 


Engine 


Aircraft Transfer Matrix 


Tc 



r*i 


A.i.B] 

throttle inputs equal 

rA 

B^column)' 

yj 


C | DJ 


C 

D-0 


x - [q(deg/sec) { a(deg) j v(kts) | 0(deg) | h(ft)| 

f 

y -[n p uo,* I n fag\ i q I a 1 V I 0 I h i y(deg)J 

» " [Zoutbd lcf< lbS ) | Z i D bd left( lbs ) | Z.nbd righ^ bs ) | Z outbd rigbfl bs )] 

u x - z(lbs) [used when all throttles have same command ] 


Lateral Dynamics 


X' 


rA | Bi 


[ X 1 


rA j b 2 I 

' X ' 



"T-- 


_ 

— 

— I- — — — — - 

— — 

yj 


[ c 1 D J 


U J 


C* D = 0 

U 2. 


where 


P./C- 


.for four engine inputs, u 


B, 


= 6 


.for one total engine input, u 2 


x ~[p(deg/sec) | r(deg/sec) j 0(deg) J </>( deg)] 


y - [A Yc ^ | p(deg/ sec) |r(deg/sec) 10 (deg) |^(deg)J 

/ 

“ [ Z outbdleftO bs ) I Z inbdleftO bs ) 1 Z inbd right 0 bs ) ! Z oulbd right 0 bs )] 

u 2 - [z(lbs)] , where z = z^,* + z^,* + (-z jnbdrigbt ) + (-z^^) 


The B matrix has four columns, each column is to be multiplied by the thrust 
input from each engine that is given in matrix u. If symmetric throttle is 
given (assume all four throttles are given the same command), the B matrix 
in longitudinal dynamics becomes a single column. Each row value in this 
column matrix Bi is equal to the sum of the corresponding row elements in 




the full order B matrix representing four engines. If differential throttle is 
given (i.e., the left engines and right engines are given same amount of 
command but in opposite directions), the B matrix in lateral dynamics 
becomes another single column matrix, B 2 . Each row element inB 2 is the sum 
of the positive value of columns 1 and 2, and the negative value of columns 3 
and 4 of each row in B. The open-loop configuration then becomes 
P=Pa/c*Peng, where Peng is the quadruple form of the engine transfer function. 

The quadruples for four different configurations were obtained as 
described in reference. 


Flight conditions for each configuration are summarized in the following 
table. 

Configuration Summary 


Config. 

Number 

Weight 

(lbs) 

BnflgEfiM 


wismkm 


1 

140,000 

4,000 

160 

0 

U P 

2 

140,000 

4,000 

145 

30 

HE 

3 

160,000 

4,000 

175 

0 

U P 

4 

140,000 

4,000 

155 

30 

U P 


The transfer functions were obtained from the quadruples using 
System Technology's Program CC. These aircraft transfer functions are listed 
here with each respective row of numbers designating the corresponding 
configuration transfer function values. The nominal configuration, number 
1, is represented by values in each row 1 below. 

Longitudinal Transfer Functions 

wiq(deg/*ec) iyq(deg/«ec) / * 

•^(ibs) l ’z(lb*) ' “long 


N r(de8) - N r(deg) / A 

”z(!b«) ™z(]be) ' “long 


igq(deg/*ec) 

'^Z(lbs) 


2.36E-04 

(0) 

(-1.17E-05) 

(0.40) 

(0.61) 

2.33E-04 

(0) 

(L4E-06) 

(0.635 , 

0.563) 

L976E-04 

(0) 

(0.292) 

(0.644) 


1.955E-04 

'W 

(2.68E-06) 

(0.819 , 

0.508) 




























iur(<kg) 

^Z(lbs) 


2.796E-05 

(0) 

(0.203) 

(0.370 , 

3.008) 

-1.819E-05 

(0) 

(0.364) 

(2.255) 

(-4.452) 

2.130E-05 

(0.167) 

(0.351 , 

a 038) 


L470E-05 

(0) 

(0.261) 

(0.460 , 1426) 



(L438E-05) 

(3.918E-02 , 

0.130) 

(0.652 , 

L382) 

(1.101E-05) 

(7.423E-02 , 

0.147) 

(0.596 , 

L375) 

(3.949E-02 , 

0.118) 

(0.649 , 

1.301) 


(L878E-05) 

(7.190E-02 , 

0.138) 

(0.588 , 

1.279) 


Lateral Transfer Functions 


|U/ 3 (deg) IU/ 5 (deg) / A 

'^Z(lb«) '^Z(lbt) ' ^lat » 


nSS > -n«*5 ) /a i 


lat ’ 


N^(deg) 

™Z(lbs) 


-L58E-03 

(-.0805) 

(.927) 

-L59E-03 

(-.0922) 

(.904) 

-L43E-03 

(-.0723) 

(.981) 

-1.44E-03 

(-.0879) 

(.940) 


U*(<kg) 

•^Z(lbs) 


3.19E-04 

(.468 , 3.65) 

2.15E-04 

(.611 , 4.17) 

2.89E-04 

(.447 , 3.96) 

Z04E-04 

(.593 , 4.33) 


(.0001) 

(1*01) 

(.116 , L05) 

(.006) 

(1.05) 

(.067 , 0.93) 

(.0028) 

(1.06) 

(.114, L08) 

(.0065) 

(1*09) 

(.060, .944) 


The engine transfer function for all configurations is given in short form 
notation by: 


G 


X(fcg) 


275 


(0.55X5) 



■ation Storage Table for Quadruples 


*Pxxxx.4U : Quadruple with four engine inputs 
*Pxxxx.lU : Quadruple with one total engine input 



































APPENDIX B: PAPERS PRODUCED IN SUPPORT OF THIS GRANT 


1. Biezad, D.J. and C.P. Azzano. "Designing Low Bandwidth Propulsive-Only Flight 
Controllers." AIAA Guidance, Navigation, and Control Conference Paper No. 91-2628CP. 
New Orleans, LA. August 12-14, 1991. 267-275. 

2. Biezad, D.J. "The Propulsive-Only Flight Control Problem." National Aerospace 
Electronics Conference, Vol. 2. Dayton, OH. May 20-24, 1991. 494-500. 

3. Biezad, D.J. and H.L. Chou. "Pilot-in-the-Loop Analysis of Propulsive-Only Flight Control 
Systems." NAECON, Vol. 2. Dayton, OH. 482-488. 

4. Biezad, D.J. and HL. Chou. "Application of QFT to the Problem of Failed In-Flight 
Controllers During Approach and Landing of a B-720 Aircraft." To be presented at 
NAECON 93. Dayton, OH. May 24-28, 1993. 


m 



DESIGNING LOW BANDWIDTH PROPULSIVE-ONLY 
FLIGHT CONTROLLERS 


Daniel J. Biezad* 

Cal Poly Slate University 
San Luis Obispo, California 

Christopher P. Azzano 
Lieutenant, United States Air Force 


Abstract 

Results from an investigation of using 
engine commands to control flight attitude are 
described. In-flight operation with simulated 
failed flight controls is reviewed and ground 
simulations of piloted propulsive-only control to 
touchdown are analyzed. A design of an optimal 
control law to assist the pilot is presented. 
Recommendations are made for more robust 
design and implementation. Results to date 
indicate that simple and effective augmented 
control can be achieved in a wide variety of failed 
configurations. 

Nomenclature 

a perturbed angle of attack (deg) 

(5 perturbed sideslip (deg) 

y perturbed flight path angle (deg) 

<p perturbed bank angle (deg) 

Tcrnd glide slope commanded (deg) 
e glide path deviation angle (deg) 

A lateral path deviation (deg) 

e T perturbed throttle (%) 

d deviation above glide path (ft) 

h altitude change - -down (ft) 
p roll rate (deg / sec) 

q pitch rate (deg / sec) 

r yaw rate (deg / sec) 

K* feedback gain for x 

G™ transfer function (s) 

Introduction 

Propulsive controls which assist 
conventional control surfaces in the attitude 
control of aircraft have been recognized as 
important enhancements of combat aircraft 
maneuverability^. In commercial operations 
such maneuverability is seldom required, but in 
the event of hydraulic failure of controls or 
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damage to control surfaces, the engines of a large 
commercial aircraft are usually capable of 
attitude control. Recent flight control failures on 
commercial aircraft, although extremely rare, 
have shown that piloted aircraft can remain 
controllable in-flight by the skillful application 
of thrust.2.3 The extreme difficulty of this task, 
however, combined with pilot stress, cannot be 
expected to result in a successful landing. 

An investigation of propulsive-only flight 
control by NASA Dryden 4 - 5 has shown it to be 
feasible for a wide variety of aircraft types and 
failure configurations. The list of aircraft flown 
include the Lear 24, Cessna 152, Piper PA-30, 
and the F-15 (single-engine aircraft required 
that the rudder be used in addition to the 
throttle). None of the in-flight tests were flown 
to touchdown. Pilot ratings were categorized by 
controlled axis and by task. Typically, 
longitudinal axis control was rated Level 2 for 
the approach and Level 3 for runway landing. 
Lateral axis control was rated Level 2 for both 
approach and landing. The pilot learning curve 
in all cases was rapid. 

Although controlled flight was always 
possible, pilots could not safely and predictably 
maneuver with the throttles alone. There may be 
sufficient control power available (presuming 
the throttles are advanced from idle), but the 
typically long time constants and couplings 
between dynamic modes make piloted flight 
precarious for demanding tasks such as landing. 
Training may alleviate the gross misapplication 
of throttles but will not guarantee safe landings. 

A pilot-assist mode which automatically 
moves the throttles in order to control attitude is 
a potential solution. Such a mode would be 
activated by the pilot in the event of complete or 
partial failure of the high bandwidth pitch and 
roll controls. This presumes that the engine 
power settings and aircraft geometry provide 
controllability under a variety of aircraft 



configurations, failure modes, and power 
settings. 

Because of the long time constants of the 
engines relative to those of the contol surface 
actuators, low-bandwidth control will be most 
effective for the long-period dynamic modes of 
the aircraft. This implies the basic airframe 
with failed controls should exhibit minimal 
stability handling qualities in flight® * 

This paper will concentrate on the major 
considerations in designing a propulsion-only 
flight control system (POFCS). The empirical 
results of ground simulations using a Boeing 720 
will be reviewed, 7 and finally an optimal linear 
design of the POFCS will be presented.® The 
paper concludes with some recommendations for 
future work. 

Ground Simulation 

Fixed-base simulations of a Boeing 720 
aircraft were performed at NASA Dryden to in- 
vestigate throttles-only control. The Boeing 
720 represents a four-engine passenger jet 
aircraft as shown in Figure 1. Asymetric thrust 
is available for roll control, but the aircraft has 
slow responding engines. Pitch-control was 
obtained by simultaneously advancing or 
retarding the throttles. A view of the simulator 
scene for approach and landing is shown in Figure 
2 . 

The Boeing 720 has a low wing with 35 
degrees of sweep. Gross attitude control in both 
the longitudinal and lateral axes during the 
simulation was possible without the use of 
electric trim. 




Low Bandwidth Control Law 

A propulsion-only flight control system 
(POFCS) must use the control power of the 
engines, assuming a stable basic airframe, to 
provide longitudinal and lateral flight path 
control under a variety of flight control failures 
throughout the flight envelope. 

Pilots must relearn how to generate lead 
compensation. There are no handling qualities 
specifications to cover this situation. The pilot 
may find it difficult to accept watching the 
throttles move with stick input The control law 
must allow pilot inputs and pilot-directed 
configuration changes without exciting large 
oscillations of the dutch roll or phugoid. 

The engine time constants must be fast 
enough to control any oscillatory mode which 
could preclude a successful landing. Relatively 
fast modes, such as the short period, must be 
stable. In other words, the configuration with 
failed controls, throughout the flight envelope, 
must be stabilizable (the uncontrollable poles 
must be stable). 9 

Boeino 720 Control Law The baseline 
configuration was gear-up, flaps-up, 10,000 ft 
pressure altitude, 160 knots, 190,000 lbs. The 
baseline control law for the four engine jet 
transport, for both the longitudinal and the 
lateral axis, was developed by trial and error in 
the flight simulator at NASA Dryden. 7 The 
baseline gains corresponding to Figures 3 and 4 
were 

{Kq.Ke.Ky.Kp.K^.Kjj} = { -4, 0,-1 ,0.5,1. 0,0.5} 


2 





As described in reference 7, ten 
configurations were then flown with the above 
set of baseline gains. The worst ratings were 
for those configurations farthest from the 
baseline weight of 190,000 lbs. Pilot 
comments for poorly rated configurations 
indicated the problem to be severe lateral 
oscillations that could not be damped 
predictably by pilot inputs. 

Classical Analysis. Linearized models of 
the longitudinal mode coupled to the glide path by 
K e are shown in Figures 5-7 and may be analyzed 
in a conventional manner as described by 
Blakelock. 10 The range to touchdown must be 
fixed for a linear analysis. Such an analysis 
shows that the baseline gains chosen are 
satisfactory longitudinally to ranges within 
1000 ft of touchdown. 

The lateral mode of coupled flight, 
however, shows an interesting feature. A two 



dimensional root locus for the lateral modes of 
response, varying Kp and K<j, , is illustrated in 
Figure 8. Note the difficulty in selecting these 
gains using conventional analysis. The lateral 
response mode has a lateral phugoid in addition to 
a dutch roll mode. Families of plots of these two 
pairs of complex roots show that varying either 
gain pushes one set of roots into the right-half 
plane. This effect of varying configurations 
exacerbates this tendency. 

Normally, given conventional flight 
controls, the pilot could compensate for this type 
of mild and slow instability. Throttles only 
control, however, even with an -augmented 
system, make such compensation extremely 
difficult for the pilot. Piloted simulations show 
that pilots are sensitive to any gain set 
significantly away from the nominal settings. 



Figure 5. Glide Slope Geometry 


Figure 6. Glide Slope Block Diagram 









pilOl 

override 



Figure 7. Longitudinal Aircraft Coupled to Glide Path 



Optimal Control Law 

A Linear Quadratic Regulator (LQR) was 
developed with modal weights chosen to provide 
robust behavior. 8 Although this choice results 
in a complex feedback structure with some loss of 
insight relative to successive loop closure, the 
opportunity to use all four engines independently 
was considered important in a flight control 
system with such degraded performance. In 
particular, this provided the capability to 
control pitch and velocity independently since the 
thrust lines of the outboard and inboard engines 
have unequal displacements along the z-body 
axis. To see this mathematically it is necessary 
to compare the eigenvalues of the controllability 
matrix (which is not done here). 


Modal Regulator Equations . The 
regulator consisted of the following feedback 


control law: 

x m = M' I AMx m + M“ , Bu (la) 

u = K m x m = K m M J x (1b) 

minimizing 

J = l/2j(x T M T QMx + u T Ru)dt (i c ) 

Q m = M T QM (Id) 


where the subscript m indicates modal coordinates 
with modal weights Q assigned directly to aircraft 

771 

dynamic modes such as the phugoid. The state 
variables for the regulator design are 

x T =[q,a,u,e,h,p,q,r,p,(p ] T (-| e ) 

where all units are radians, feet, and seconds, 
and where the control is given by 

u = [four throttles] (If) 

The LQR design condition was 4,000 ft MSL, 
175 KCAS, 160,000 lbs, with gear and flaps up. The 
plant, throttle control, and feedback matrices for this 
condition (eg 20.85% MAC) are given in Figure 9. 
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Figure 9. LQR Plant, Control, and Feedback Matrices 
for Design Configuration 


The insight gained by using modal cost 
weights can effectively be seen by comparing the 
open and closed-loop responses of sideslip to an 
Initial sideslip of 10 degrees, as shown in Figure 


10. One of the penalties of this approach, 
however, was the normalization of units so that 
weights in the cost function did not differ by 
many orders of magnitude. 




Figure 10. Cost Function Weights Trade-offs (Spiral/Dutch Roll) 
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Increasing weights on the dynamic 
modes provided tighter control, but performance 
suffered when the configurations were altered. 
This issue of robustness did not exist for the 
short period or roll modes since they were 
difficult to realistically excite by the engines, 
especially when engine lag was taken into 
account. 

The final gain matrix in Figure 9 
eliminated state feedbacks for altitude, h, and for 
the engine model states described in the section 
below with no adverse effects. 


Engine Model. The values in the B 
matrix were obtained from steady-state 
perturbations in response to thrust. In order to 
use % throttle position, ej, and not pounds of 
thrust in the control u of Equation (If), the A 
matrix was augmented with a second order engine 
model of the form 


T 

f 


' 0 1 ■ 

rr 


o 1 

-co 2 . -2 £co n _ 

T 

L J 

+ 

i 

1 


( 2 ) 


where co 2 . =2.5, £=.802, and k t„=250. These 
parameters were chosen to match engine 
transients at a nominal steady-state throttle 
setting of 20%. The transfer function form of 
the engine transient was 

g t 100- 0.55-5.0 
* (s + 0.55 )(j + 5.0) 

The conservative engine model of 
Equation{2) overestimated the gain K Tm so the 


controller would not de-stabilize the system by 
acting on a low authority plant. Equation (3), on 
the other hand, is accurate for most steady-state 
conditions on approach. 

System Dynamics. The aircraft and 
engine models were normalized as described 
above by the factors shown below in Table 1. 

Table 1. Normalization Factors 


Dimension 

Units 

Factor 

Angle 

Radians 

.001 

Force 

Pounds 

5906 

Distance 

Feet 

1 

Throttle 

Percent 

0.1 


The weights for the dynamic modes given the 
normalization factors are shown in Table 2. 


Table 2. LQR Weights 


Mode 

Weight 

Period (sec) 

£ 

Short Period 

1 

10 

0.7 

Phugoid 

10 

57 

0.01 

Dutch Roll 

200 

6.1 

0.10 

Spiral 

0.5 

68 

... 

Roll 

1 

1.2 

... 


The discrete version 1 1 of this system 
at 50 Hz has the open and closed-loop system 
roots as shown in Figure 11. 


0 









* r * * * , 

.. 

. • - . \ 

1 . :-_c> 



-0 .2 




y///; 

- 


-0.4 

1*. . . , 

\ 

d . , 9 .pon — }<?<? p • ■ r b 6 1 s . . ** 

.•** X*' y 

’ t 

.v / 





x 'ploBdd-lijop robts 

.♦* / 
,'k. 




“0 .6 


-J — — - i i - * * ' 





0 

.2 0.4 0.6 

0.8 

1 


Figure 11. 

Discrete System Roots 

at 50 Hz 




6 






Pilot Command Interface. A command 
interface was designed into the LQR loop based on 
a pseudo rate command for both pitch and roll. 
This allowed a "batch" test of the linear system 
prior to implementing the control law on the hi- 
fidelity nonlinear simulation. 

Control authority that would hot hinder or 
"wash out” pilot commands was provided by 
translating the command into a pseudo rate 
command. This rate command was digitally 
integrated over time to determine pitch and bank 
attitude command. Limiters were also inserted to 
prevent saturation. Parameter values are shown 
in Table 3 for the final LQR controller, including 
pilot interface, of Figure 12. - 


Table 3. Pilot Interface 


Parameter 

Value 

Qstick 

pitch limit 

T s 

Psiick 

bank limit 

0.25 

1 6, 000/ gross wt 
0.02 
-0.25 
6 1 0 6 / 1 22 


Note that the lateral commands take 
precedence over longitudinal ones, emphasizing 
the rationale that survivability depends 
primariy on wings-level flight and touchdown. 
Also, setting limits as a function of throttle lever 
position will likely be an impractical 
implementation. The independent control of 
velocity and pitch attitude was not accomplished 
in this investigation. 

In general, "batch" linear simulations 
prior to piloted simulation predicted higher gains 
for adequate control than were required. 

Unstable pilot-in-the-loop operation in the 
nonlinear simulation required that the weights be 
reduced to those shown in Figure 9. Those gains 
provided adequate closed-loop performance on the 
piloted simulator. 

Linear system performance did not model 
the coupled longitudinal and lateral modes, engine 
nonlinearities, and unequal engine spool-up and 
spool-down times. The nonlinear simulation also 
exhibited more Dutch Roll damping. Piloted 
simulation results follow in the next section. 
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Figure 12. LQR Propulsion-Only Control with Pilot Interface 
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Conclusions 


Piloted Evaluation. The LQR 
propulsive-only controller of Figure 12 was 
implemented on a high-fidelity nonlinear 
simulator at NASA Dryden and evaluated by test 
pilots. Stick deflection produced thrust 
commands that were observed on the engine 
instruments but which did not physically move 
the throttles. 

Pilots were asked to evaluate changes in 
altitude, velocity, heading, ground tracK, and 
flare performance. The LQR implementation did 
not allow throttles to be physically moved to 
control velocity since the LQR controller would 
interpret this as a pitch command. Such motion 
would also alter the bounds in the limiters shown 
in Figure 12. 

Up-and-away maneuvering did not 
require excessive pilot workload. The 
longitudinal implementation, however, made it 
difficult to fly level, especially when rolling out 
of turns. Because of the many development 
changes which occurred during piloted flight, 
pilot ratings were not documented. In general the 
aircraft could be considered Level 2. These 
ratings are similar to those described in this 
paper's section on the Boeing 720 control law. 

The qualitative results of the piloted 
evaluations are summarized in Table 4. As 
ex pec ted for this type of controller, performance 
degraded when the LQR control law was 
implemented on different failure configurations. 


Table 4. Qualitative Pilot Evaluation 


Task 

Pilot 

Comment 

Altitude 

Change 

'^eFocity” - 

Change 

Holds a rate of climb, but 
return to level is difficult. 

Somewhat "mysterious" 

Heading 

Change 

"FfoTd 

Sink Rate 

■fiords sfeady - staTelurn welC 
but roll -rate command more intuitive 
than bank angle command Difficult 
to maintain altitude. 

"Accepfable^after - learriihg'curve" "Flight " 
path angle command is preferable 

"Flo7d 

Ground 

Track 

Acceptable, but lightly damped roll 
is somewhat bothersome. 

Flare 

Too mucHTag! Wilfrequire some 
practice to determine when to 
initiate the flare. 


Simulations of a Boeing 720 aircraft with 
failed flight controls show that a propulsion-only 
flight control system (POFCS) is feasible. 
Classical analysis using successive loop closure 
results in simple, effective controllers. Two 
lightly damped lateral modes, however, can 
become unstable given minor gain variations or 
changes in configuration. This suggests that 
compensation should augment the gains to provide 
a more robust and stable system. 

A LQR augmentation scheme designed using 
optimal control was flown successfully under 
pilot control but was not a significant 
improvement over the gains set by classical 
analysis. The design implementation employed 
pseudo rate commands, required limiters, and did 
not allow the pilot to use the throttles for 
velocity control independent of pitch. 

Despite these limitations it was 
demonstrated by piloted evaluations that the 
POFCS concept is feasible and may be 
implemented as a back-up pilot assist mode when 
normal flight control has failed. If an optimal 
controller is employed, an improved pilot 
interface will be required as well as provisions 
for velocity control independent of pitch attitude. 
In particular, the use of differential inboard- 
outboard thrust should be investigated as a way to 
uncouple the velocity and pitch modes. 

Pilots indicated that improved handling 
qualities are desired. The wandering bank angle 
and difficulty in achieving level flight should be 
rectified. Pitch sensitivity may be reduceo by 
transition to a flight path angle command system. 

Research is in progress to investigate the 
use of compensators for failed flight control 
configurations. Better performance and more 
robust behavior for off-design failure conditions 
are desired. Flight operations which are coupled 
to the glide path will be analyzed, and the 
resulting controller is intended to be tested in 
both ground and inflight simulations. 
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The Propulsive-Only Flight Control Problem 

by 

Daniel J. Biezad 

Cal Poly State University 
San Luis Obispo, CA 93407 


Abstract 

Altitude control of aircraft using only the 
throttles is investigated The long time constants of 
both the engines and of the aircraft dynamics, 
together with the coupling between longitudinal and 
lateral aircraft modes, make piloted flight with failed 
control surfaces hazardous, especially when 
attempting to land. This research documents the 
results of in-flight operation using simulated failed 
flight controls and ground simulations of piloted 
propulsive-only control to touchdown. Augmentation 
control laws to assist the pilot are described using 
both optimal contrc' and classical feedback methods. 
Piloted simulation using augmentation shows that 
simple and effective augmented control can be 
achieved in a wide variety of failed configurations. 

Nomenclature 

q perturbed pitch rate (rad/sec) 

a perturbed angle of attack (rad) 

u perturbed velocity (ft/sec) 

6 perturbed pitch angle (rad) 

h altitude change-down (ft) 

y perturbed flight path angle (rad) 

T glide slope deviation angle (rad) 

Kq.Kg.Ky.Kr long, feedback gains 
p perturbed roll rate (rad/sec) 

r perturbed yaw rate (rad/sec) 

P perturbed sideslip (rad) 

<t> perturbed bank angle (rad) 

o lateral offset angle from runway 

K<t>,Kp.Kp,Ko lateral feedback gains 
Syp perturbed throttle {%) 

Syp perturbed thrust 

5 D y perturbed differential thrust 

Introduction 

The failure of hydraulic power to primary flight 
control systems is an extremely rare occurrence in- 
flight. Such failures have occurred, however, and their 
consequences have been especially tragic in 
commercial operations. In a few of these failures the 
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aircraft remained controllable in-flight by the skillful 
application of thrust, but the extreme difficulty of this 
task, combined with the stress of the emergency, did 
not allow a successful landing. 

The fundamental problem is that an aircraft 
cannot be easily and predictably maneuvered by the 
pilot with the throttles alone. Although the control 
power is often sufficient to fly the aircraft, the long time 
constants and couplings between dynamic modes 
make pilot control uncertain and precarious for 
demanding tasks such as landing. Exposure to these 
situations in training simulations may alleviate the 
gross misapplication of throttles, but will not eliminate 
the potential for a serious accident to occur. 

A complimentary solution to more training 
would be the addition of a simple, low cost pilot-assist 
mode to be activated by the pilot in the event of 
complete failure of the high bandwidth pitch and roll 
controls. The goal of such a system would be to 
provide acceptable flying qualities by driving the 
throttles through pilot command inputs from the 
control column. Although the requirements are many 
for such a system, three factors are worthy of note. 

First, the engine power settings and 
mounting geometry must provide controllability in a 
mathematical sense for the aircraft equations of 
motion under a variety of aircraft configurations and 
failure modes. Second, and possibly the most 
important factor, the low-bandwidth control system will 
be coupled with a pilot who will be stressed and 
anxious under actual emergency conditions. There 
are no handling qualities requirements to guide the 
designer here, and even those in military 
specifications 1 are inapplicable. Third, because of the 
long time constants of the engines relative to those of 
the control surface actuators, low-bandwidth control 
will be most effective for the long-period dynamic 
modes of the aircraft. This poses a special difficulty 
for stabilizing the lateral Dutch Roll mode, which may 
not be a slow mode relative to the engine response. 

This paper will concentrate on the 
development of a practical propulsion-only flight 
control system (POFCS). The first major issue of 
controllability due to engine power and geometry will 
be broadly surveyed in the next section titled “In- 
Flight Simulations". 2.3 Next, the empirical results of 


ORIGINAL. PAGE fS 
OF POOR QUALITY 




ground simulations using a Boeing 720 aircraft model 
with failed flight controls will be presented. The third 
section will present development issues for the 
POFCS and highlight the difficulties of achieving a 
robust, practical design. 


In-Flight Simulations 

Preliminary investigations of throttle-only 
aircraft control in-flight have been conducted by 
NASA Dryden Flight Research Facility. The list of 
aircraft flown includes the Lear 24, Cessna 152, Piper 
PA-30, and the F-15. Single-engine aircraft required 
that the rudder be used in addition to the throttle. 

None of the in-flight tests were flown to touchdown. 
Pilot ratings could be .categorized by controlled axis 
and by task. Typically, longitudinal axis control was 
Level 2 for approach and Level 3 for landing on a 
runway. Lateral axis control was Level 2 for both 
approach and landing. The pilot learning curve in all 
cases was rapid. 

F-15throt1les-onlv. The basic F-15, shown in 
Figure 1 , has a high wing with approximately 45 
degrees of sweep and two vertical tails. Although the 
engines are near the aircraft centerline, flight tests 
showed roll rates from differential throttles up to 5 
deg/sec over a significant part of the flight envelope. 
Pitch control using throttles alone was available but 
inadequate below 200 knots. 

With yaw augmentation systems off, three test 
pitots found that differential thrust atone provided 
good bank angle control as well as roll rate response. 
With nose up trim, and stick and rudder centered, they 
were able to exercise crude altitude and heading 
control. 


moment from a thrust increase. The aircraft tnus n au 
to be flown using its inherent speed stability, leaving 
an undamped phugoid and a pitch rate response of 
approximately 0.2 deg/sec. 

Test pitots tried to use bank angle to damp the 
phugoid with inconsistent results. Only when electric 
pitch trim was made available could a tractable 
approach be flown. As with all in-flight tests, no 
landings were attempted using throttles atone. 



Cessna 152. The single-engine light trainer 
shown in Figure 3 has a high wing and conventional 
tail. Rudder was required for directional control, but 
the throttle provided adequate control of the steady- 
state speed stability of the aircraft. 

Phugoid excitation required damping by pitot 
application of throttle which was unnatural but easily 
learned. Pitots stated that they could have landed the 
aircraft using throttle and rudder alone. 



Lear 24 throttles-onlv. The twin engine 
executive jet shown in Figure 2 has a T-tail and fast 
responding engines. In-flight thrust conlrol of roll rate 
was effective, reaching 20 deg/sec near 250 knots. 
Pitch response due to thrust was very poor, due to 
the high engine location which caused a nose down 



Figure 3. Cessna 152 
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Piper PA-30. This twin-engine aircraft shown 
in Figure 4 has a low wing and conventional tail. The 
roll control power was considerable but very non- 
linear, requiring extensive pilot adaptation. Roll rates 
were observed near 10 deg/sec but bank angle 
control was very difficult. 

Pitch control from throttles-alone came from 
the inherent speed stability of the aircraft. Pilot 
damping of the phugoid was difficult and would have 
made landing under throttles only control dangerous. 
Providing electric pitch trim alleviated the problem to a 
great degree, and it was possible for two pilots to 
simulaneously control flight in this manner. 



Figure 5. Boeing 727 


Ground Simulations 

Full six degree-of-freedom simulations of 
large aircraft were performed at NASA Dryden to 
investigate throttles-only control. The Boeing 727 
and the Boeing 720 represent three- and four-engine 
variants of passenger jet aircraft as shown in Figures 5 
and 6. Although considerable differential thrust exists 
for roll control, both aircraft have slow responding 
engines making damping of the dutch roll and 
phugoid modes difficult. Differential thrust was not 
used to control pitch attitude. A view of the simulator 
scene for approach and landing is shown in Figure 7. 



Boeing 727. This three engine transport has 
a swept engine and a T-tail. From a level flight trimmed 
condition, throttles moved in concert produced about 
0.5 deg/sec in pitch, and throttles moved differentially 
produced 3 deg/sec roll response. Electric trim was 
required to damp the phugoid sufficiently for a landing 
on a "field." Without the trim extensive practice was 
required (over 2 hours). 

When two pilots divided the control task by 
axis, they could successfully land the aircraft but not 
on a runway. Considerable care was required not to 
excite the dutch roll and phugoid modes. Pilots found 
this unnatural and especially difficult when 
approaching touchdown. Level 3 ratings require 
some sort of stability augmentation for safe flight. 

Boeinn 720. This four engine jet transport 
has a low wing with 35 degrees of sweep. Gross 
attitude control in both the longitudinal and lateral 
axes was possible without the use of electric trim. If 
plots split the tasks and used electric trim, a runway 
landing could be made. 



Figure 7. Simulation Visual 


An augmentation system to control flight path 
angle, y, was developed to convert conventional pilot 
control inputs into throttle motions. Using this system 
a single pilot could successfully land the aircraft with 
practice on a simulated runway. There still was a 
tendency for pitot-induced oscillations (PIO) near the 
ground, especially in the lateral axis, and control 
inputs for stability were requied to be very small to 
avoid exciting the oscillatory modes. 
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Control Law Development 

An ideal propulsion-only flight control system 
(POFCS) would provide acceptable handling qualities 
in the event of any type of flight control malfunction. 
Essentially, however, the implemented control law 
must perform three functions. First, it must use the 
control power of the engines, in concert with the 
stability derivatives of the aircraft, to provide 
longitudinal and lateral flight path control under a 
variety of flight control failures throughout the flight 
envelope. This is difficult to achieve when the aircraft 
is far from a trimmed condition, or when the aircraft 
must descend rapidly under conditions of low thrust. 

Second, the control law must allow pilot inputs 
and pilot-directed configuration changes without 
exciting large oscillations of the dutch roll or phugoid. 
Pilots who are trying to perform a high gain task with a 
low bandwidth control , such as landing an aircraft, 
must relearn how to generate lead compensation. 
There are no handling qualities specifications to cover 
this situation. The pilot must accept watching the 
throttles move with stick input and must resist the 
natural tendency to pull power off during the flare. 
Such a control law must be integrated with actual pilot 
inputs at each step during its development. 

Finally, the engine time constants must be 
fast enough to control any oscillatory mode which 
could preclude a successful landing. Other relatively 
fast modes, such as the short period, must be stable. 
This condition simply states that the configuration with 


failed controls, at points in the flight envelope from 
landing to cruise conditions, must be stabilizable (the 
uncontrollable poles must be stable). 4 

Boeing 720 Control Law The baseline 
control law for the four engine jet transport, for both 
the longitudinal and the lateral axis, was developed by 
trial and error in the flight simulator at NASA Dryden. 
The augmentation control law for each axis was 
developed from a baseline aircraft configuration (gear- 
up, flaps-up, 10,000 ft pressure altitude, 160 knots, 
190,000 lbs). For this baseline Level 2 pilot ratings 
were recorded for each axis given the task of landing 
on a runway. The baseline gains were 

{K q ,K0,Ky,Kp,K<j > ,Kp} = { -4.0.-1.0.5.0.5.1.0} 

where other potential gains, such as K r , were tried but 
not kept in the baseline since the pilots noted no 
significant improvement in aircraft response. 

Ten configurations were then flown with the 
above set of baseline gains and the results are 
summarized in Table 1. The worst ratings, those for 
configurations 5, 8, and 9, were for those 
configurations farthest from the baseline weight of 
190,000 lbs. The poor rating (8) for configuration 3 
can be considered an anomaly based on pilot 
comments of full stick throw and starting in too close. 
Otherwise, pilot comments for poorly rated 
configurations indicate the problem to be severe 
lateral oscillations that cannot be damped predictably 
by pilot inputs. 


Table 1. Cooper-Harper Ratings for Boeing 720 Simulation to Touchdown 
Baseline: 190000 lbs, 10000 ft, gear-up, flaps-up, 160 knots, fight turbulence 
{Kq.Ke.Ky.Kp.K^Kp} = { -4.0.-1.0.5.0.5.1.0} 


CONFIGURATION Pilot Comments 

# Altitude Weight Speed Gear/ Rating 
(ft x 1 0 3 ) (tb x 1 0 3 ) (knots) Raps 


1 . 

4 

140 

160 

UpUp 

5-6 

Landed long; needed small input 

2. 

4 

190 

160 

UfVUp 

4-5 

OK! 

3. 

10 

190 

160 

Up/Up 

8 

In too tight; kept hitting full stick thrdw; no control power 
ooukJ not get back from 30° bank; excited phugoid. 

4. 

10 

140 

160 

UpUp 

6 

Kept VVI above 500 fpm to keep control power; OK! 

5. 

4 

140 

160 

Up60° 

9.7 

Crashed! Initial excitation caused undampabte dutch roll. 
Second try better but could not get in loop safely. 

6. 

4 

190 

160 

UpV50° 

7 

Could not get in control loop safely laterally; pulsed inputs. 

7. 

4 

190 

130 

Up60° 

5 

Small inputs required; excellent control; could damp roll. 

8. 

4 

140 

130 

Up50° 

8 

Controlled crash off of runway!; open loop only; very difficult 
to damp dutch roll; pulsed inputs did not hefc>. 

9. 

4 

140 

130 

Down/500 7 

Controlled touchdown off runway; same comments as #8. 

10 . 

4 

190 

130 

D0W1V5QO 5 

Excellent! Same comments as #7. 


4 




Classical Analysis. A linearized model of 
configuration 1 of Table 1 is shown in Figure 8 . A 
two dimensional root locus for the lateral modes of 
response, varying Kp and K<f> , is illustrated in Figure 9. 
Note the difficulty in selecting these gains using 
conventional analysis. The lateral response mode has 
a lateral phugoid in addition to a dutch roll mode. 
Families of plots of these two pairs of complex roots 
show that varying either gain pushes one set of roots 


Engine Representation 


into the right-half plane. This effect of varying 
configurations exacerbates this tendency. 

Normally, given conventional flight controls, 
the pilot could compensate for this type of mild and 
slow instability. Throttles-only control, however, even 
with an augmented system, make such compensation 
extremely difficult for the pilot 
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Figure 8. Linearized Model of Configuration 1 
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Optimal Control Analysis. An optimal control 
law was developed by Azzano for the linearized 
system shown in Figure 10. In his development a 
pilot command is directly translated into a single state- 
variable command. All four throttles were 
independently controlled, taking advantage of the 
vertical offset of the engines to allow pitching 
moments independent of airspeed change. 

The design condition for the optimal controller 
was 4000 ft altitude, 160000 lbs, 175 knots, gear-up, 
and flaps-up. The engine lag was modeled as a 
second order system. Azzano’s final design, 
incorporating many limiters, is shown in Figure II, 
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Figure 9. Lateral Root Locus 
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Figure 10. Matrices for LQG Model 
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Figure 11. Linearized LQG Model for configuration 1 at 175 knots 


Although open- and closed-loop responses 
of the ground simulator and the linear system matched 
reasonably well for small inputs, pilots found the 
optimally-designed augmentation system difficult to 
use in a few key areas. Pilots described returning to 
level flight coming out of a turn as unpredictable and 


velocity change as ’mysterious'. The flare maneuver 
was not responsive enough. They typically 
complained of a "wandering bank angle" and had a 
tendency to excite pilot induced oscillations even 
when the controller gains were reduced by changing 
the weights on the cost function. 
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The augmentation system designed using 
optimal control theory performed poorly when other 
failed configurations were tested. This was expected 
for this type of controller. 


Conclusions 

The results of in-flight operation using 
simulated failed flight controls for a variety of different 
aircraft show that the throttles can be effective low 
bandwidth controllers. When throttle controllability is a 
problem, electric pitch trim and/or rudder input may be 
required for safe flight. Pilot learning is rapid, but 
performance during high gain tasks such as 
touchdown is not predictable and requires 
augmentation. The pilot has difficulty generating lead 
compensation for a low bandwidth, lightly damped 
control system. 

Ground simulation of a Boeing 720 four- 
engine jet transport showed that a simple 
augmentation system could assist the pilot. Problems 
were apparent, however, when the failed 
configuration being flown deviated from the design 
failed configuration, with pilot rating being most 
sensitive to aircraft weight and center of gravity. 

Classical analysis highlighted the problem to 
be two lightly damped lateral modes which became 
unstable if either the roll rate or the beta feedback 
gains were increased from a nominal setting. Flight off 
of the design condition exacerbated the problem. 

The pilot found it difficult to control the off-design 
flight configuration in the lateral mode. 

An augmentation scheme designed using 
optimal control was successful under pilot control but 
required many limiters and adjustments of design 
weights.’ The performance away from design 
condition was not acceptable. 

Work is in progress to use quantitative 
feedback theory to build a simple compensator for the 
failed flight control configurations which will be more 
robust with respect to off-design conditions. 

Acknowledgment 

All of this work was made possible by the 
support of technical personnel at NASA Dryden, 
especially Bill Burcham, Glenn Gilyard, and Joe 
Conley. Use of the simulation facility was vital to the 
completion of this investigation. 


References 

1. MIL-STD-1797A, “Flying Qualities of Piloted 
Vehicles", limited distribution, ASD/ENES, Wright- 
Patterson AFB, Ohio, 30 Jan 1990. 

2. Burcham, F., Fullerton, G., Gilyard, G., Conley, J. 
and J. Stewart, "A Preliminary Look at Propulsion- 
Enhanced Flight Controls Research," NASA Dryden 
internal document, 1990. 

3. Gilyard, G., Conley, J., Le, J., and W. Burcham, "A 
Simulation Evaluation of a Four-Engine Jet Transport 
Using Engine Thrust Modulation for Flight Path 
Control," 27th Joint Propulsion Conference (to be 
presented), June 24-26, 1991, Sacramento, CA. 

4. Kwakemaak, H. and R. Sivan, Linear Optimal 
Control Systems. John Wiley and Sons, 1^72. 

5. Azzano, C.P., "A Preliminary Look at Optimal Multi- 

Variable Design of Propulsion-Only Flight Controllers 
for Jet Transport Aircraft," NASA Dryden internal 
document. Sept 21, 1990. .. 5 

6. Roskam. J.. Airplane Flight Dynamics and 
Automatic Flight Controls. Roskam Aviation and 
Engineering Corp., Ottawa, Kansas, 1979. 

7. McRuer, D., Ashkenas, I., and D. Graham, Aircraft 

Dynamics and Automatic Control. Princeton University 
Press. 1973. . - " , -• 


7 


Pilot-in-the-Loop Analysis of 
Propulsive-Only Flight Control Systems 

b y 

Hwei-Lan Chou 
Daniel J. Biezad 

Cal Poly State University 
San Luis Obispo, CA 93407 


ABSTRACT 

Lomitudinal control system architec- 
tures are presented which directly couple 
flight stick motions to throttle commands for a 
multi-engine aircraft. This coupling enables 
positive attitude control with complete failure 
of the flight control system. The architectures 
chosen vary from simple feedback gains to 
classical lead-lag compensators with and 
without prefilters. Each architecture is 
reviewed for its appropriateness for piloted 
flight. The control systems are then analyzed 
with pilot-in-the-loop metrics related to 
bandwidth required for landing. Results 
indicate that current and proposed bandwidth 
requirements should be modified for throttles 
only flight control. Pilot ratings consistently 
showed better ratings than predicted by analy- 
sis. Recommendations are made for more 
robust design and implementation. The use of 
Quantitative Feecback Threory for compensator 
design is discussed. Although simple and 
effective augmented control can be achieved in 
a wide variety of failed configurations, a few 
configuration characteristics are dominant for 


pilot-in-the-loop control. These characteris- 
tics will be tested in a simulator study involving 
failed flight controls for a multi-engine 
aircraft. 


NOTATION 

q pitch rate (deg/sec) 

a perturbed angle of attack (deg) 

u perturbed velocity (ft/sec) 

9 perturbed pitch angle (deg) 

h altitude change-down (ft) 

Y perturbed flight path angle (deg) 

r glide slope deviation angle (deg) 

roii! 

°in transfer function (s) 

(a) short form for (s+a) 

(C,ro„) short form for s’ + 2 Co s + © J 
Kq pitch rate feedback gain 

K y gamma loop feed back gain 

&rc throttle command (%) 

z thrust (lbs) 

<5 S stick input (full deflection=1 unit) 

P quadruple state-space representation 

M u dimensional speed derivative 

M a static stability dimensional derivative 


1 



INTRODUCTION 


THROTTLES-ONLY SYSTEM SURVEYS 


Work at NASA Dryden has shown that 
compensated thrust modulation coupled to 
Might stick motion provides a positive degree 
of flight controllability in the event of 
complete failure of the flight control system. 
Feedback control laws developed empirically 
had dramatically improved the pilot ratings 
from Level 3 to Level 2 for the simulated 
approach and landing of a Boeing 720 with 
failed flight controls 13 . Initial work on the 
modeling of these control systems showed that 
relatively simple feedback architectures, as 
well as those based on optimal control theory, 
could ease the piloting task for throttles-only 
flight unless moderate turbulence was 
encountered. 4 5 

The main thrust of research reported 
here has been to investigate the effect of 
throttles-only flight control on the flying 
qualities of multi-engine aircraft. Analytical 
system surveys are accomplished to explain 
this improvement from a handling qualities 
point of view. The pilot-in-the-loop metrics 
used in the investigation are primarily related 
to bandwidth criteria as reported in the 
literature. 6 

Previous work was extended by devel- 
oping classical compensator designs with and 
without prefiltering to further improve the 
piloted ratings. The design goal was to find a 
robust controller for throttle-only control 
under various approach and landing flight 
conditions. Designs obtained from optimal 
control theory showed performance sensitivity 
to configuration changes 5 . 

All work assumes that the aircraft 
configuration has a positive My dimensional 
derivative and positive stability (M a < 0). 
System surveys follow, then the design archi- 
tectures are analyzed. An expanded Appendix 
describes the aircraft configurations. 


The basic system model as shown in the 
Appendix has four variations of configuration. 
The engine and bare airframe state-space mod- 
els, called quadruples 7 , were derived from 
perturbations of the full non-linear equations of 
motion about trim. Transfer functions used in 
design were then approximated with low order 
fits over the frequency range of effective throt- 
tle control. 



The spool-up and spool-down engine 
dynamics for the B-720 engine are shown in 
Figure 1. The empirical transfer function 
developed is given in short form notation by 
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275 

(0.55X5) 


The above equation is illustrated in Fig. 2 
over low frequency ranges up to 1 .0 rad/sec. 
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Figure 3 Pitch Rate to Thrust Bode 


Bare Airframe. It is apparent from the 
engine bode diagrams that severe bandwidth 
attenuation occurs beyond frequencies of 
1 rad/sec. It may not be possible, therefore, 
to increase the closed-loop bandwidth beyond 
1 rad/soc within the range of available thrust 


This can be seen in the pitch rate “q” to 
thrust "z” transfer function of the bare 
airframe shown in Figure 3. The full-order 

thansfer function shows that 80 db of 

gain must be added to yield a crossover 
frequency beyond 1 rad/sec. This corresponds 
to 10,000 lbs of full thrust from each engine, 
which would not be practical for approach and 
landing. 


a i , pq(dc*/itc) 

A low order fit to b , c ^, is also depicted in 
Figure 3 and is very accurate near the phugoid 
frequency. Piloted flight of the unaugmented 
aircraft was consistently Level 3. The main 
difficulties were the lightly damped phugoid 
and the low bandwidth throttle control. The 
open-loop response of pitch angle to a full 
deflection step stick input is shown in Figure 4 
with all compensation set to unity (see 
Appendix). 



The accuracy of the low order fit near the 
phugoid frequency means that, to a first order 
approximation 9 , the phugoid frequency and 
damping are found from 


2 C®. = -X u + 


M u (X a -g) 


2 M r 


cu = 


LI 


and for conventional transport aircraft can 
shown to be roughly proportional to M u . 


be 


It should be strongly noted here for the 
classic case of M u =0 and for negative values of 
M y (Mach luck) that the aircraft cannot be 
practically flown with throttles alone unless 
rotational control in pitch is added. Difficulties 
will also be encountered as M a becomes small 
(aft eg location). Both of these cases require the 
addition of an effective rotational controller 
about the pitch axis. This may be achieved using 
differential inboard and outboard thrust, 
provided the inboard engines are a different 
distance from the aircraft xy-plane than the 
outboard engines. These configuration charac- 
teristics determine the innate capability for 
throttles-only piloted control. 
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Figure 5 Pilch Angle System Survey 

FEEDBACK ARCHITECTURES 


The generic feedback architeclure is 
given in the Appendix. An effort was made in 
the designs to keep the structure simple, and so 
in all cases the flight path compensation was 
unity. The open-loop pitch angle to stick root 
locus. Bode, and "Siggy" plots 8 are shown in 
Figure 5. They are characterized by excessive 
resonance at co np h. low phase and gain margins, 
low crossover frequency, and large phase angle 
roll-off. The open-loop (OL) system is 

G W '2 s 8.32(0.4)(0.61) 

^ (.0039, 0.1 3)(0.65, 1.38)(0. 55)(5.0) 


The root locus of the open-loop system 
makes it apparent that any feedback is limited 
by the phugoid roots going unstable. 


Empirical Feedback. This longitudinal 
control law was developed by trial and error in 
the simulator at NASA Dryden with a pilot in 
the loop. It is given by 


i nXm (d> g) 
*■ (units) 


nS u (%) 


,K q ,K r ,G 


e tn < dc 8h 
r(d c g ) f 


(10,10,4,1,1) 


The system survey for this structure is 
shown in Figure 6. 
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Figure 6 Survey with Empirical Feedback 


It can be seen that the q and y feedback 
loops removed the resonance at the phugoid 
frequency along with the rapid phase drop. The 
gain and phase margins, however, are still low. 
The q loop closure caused the increase in phugoid 
damping, and the y loop closure provided an 
additional 70% increase in settling time. The 
empirical feedback essentially cancelled the 
modified engine mode at -0.397 as shown below. 


G™ . 

o t (units) 


loops 

dosed 


8. 42(0.4X0.61) 

(0. 5 1 8, 0. 244 )(0. 5 1 7, 1.5X0. 397)(5. 1 6) 


Classical Feedback Design. Classical 
compensation was designed to address the low 
gain and phase margins and to increase system 
bandwidth within the practical limits of the 
throttle command. The compensation chosen was 


(deg) »~5 It (5>) y V n fl '" (dc£) i 
uniis)» u e,(deg)> lv -q» , 'j.» u j.(dcg) / 

i a (S + 0.55) (S + 0.65) 

(S + 0.65)’ ’ (s + 1.3) ’ 



The survey for this system is similar to Figure 
6 and not repeated here. 
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The classical design improved the 
empirical one primarily by increasing the 
phase margin of the pitch angle to stick 
transfer function from 13 to 26 degrees. The 
crossover frequency remained near 

0.98 rad/sec. and the steady state performance 
increased 1 0%. — 

The improvement in phase margin made 
the controller more robust when used to fly the 
other configurations. The empirical controller 
was also surprisingly robust when used to fly 
the other configurations. A complete discus- 
sion of this is found in Reference 9. 

Further improvements in bandwidth 
could be achieved only by substantially raising 
the compensator gain. This resulted in exces- 
- sive control (thrust). A comparison of the 
thrust response to a full stick step deflection 
for the different feedback architectures is 
given in Figure 7. It was assumed that the 
throttle command could be moved instanta- 
neously. A throttle actuator would introduce an 
additional lag. 

Compensators currently being designed 
using Quantitative Feedback Theory are having 
similar difficulty meeting reasonable limits on 


control activity when the design closed-loop 
bandwidth is near 1 rad/sec. A design procedure 
is being developed to determine the achievable 
closed-loop bandwidth for a set of configurations 
given a bandwidth limit on a primary controller. 

CONCLUSIONS 

Bandwidth requirements on pitch to stick 
response should reach 3 rad/sec for acceptable 
pilot ratings 6 . Augmented throttles-only flight 
could not reach beyond 1 rad/sec, and received 
acceptable Level 2 ratings unless moderate 
turbulence was applied to the simulation. Work 
in progress at Systems Technology Inc. is estab- 
lishing bandwidth limits for large, landing 
aircraft, and these limits will be used to design 
future compensators. Within the limits set by 
key configuration variables Mu and M u , simple 
classical compensators that increase the phase 
margin result in acceptable pilot ratings for 
throttles-only flight. 
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APPENDIX: B-720 CONFIGURATIONS 

The B-720 piloted simulation can be represented by the following block diagram: 



The ‘aircraft above represents both the engine and the bare airframe dynamics. The engine is ap- 
proximated by a transfer function and the bare airframe dynamics are represented mathematically by a 
single quadruple, Pa/c, shown as follows. 


Throttle command 
% rpm 

< 5 Tc (%) 


Engine 


Aircraft Transfer Matrix 



y 


‘A ! B' 

throttle inputs equal 

A IB, (column)' 

C \ D 


C ! D = 0 


x = [q(deg/sec) J a(deg) \ v(kts) | 0(deg) ! h(ft)] T 
y=[*Ws ! n gs ; q : « ; v j <9 ; h i r(deg)] T 
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The thrust control depends on whether four engines receive independent or identical commands: 


u [ z ombd left Obs) | z mM]cfl (lbs) ] z jnMfighl (lbs) j z.HBMrigh.Gbs)] 

u, = z(lbs) [used when all throttles have same command] 


Note when all four throttles are given the same command from the pilot stick input the B 
matrix becomes a single column. Each row value in this column matrix is equal to the sum of the 
corresponding row elements in the full order B matrix representing four engines. The open-loop 
configuration then becomes P=Pa/c*Pe, where Pe is the quadruple form of the engine transfer 

pZ(IbO 

function, <*(»>. The quadruples for four different configurations were obtained as described in Refer- 
ence 9. 


The flight conditions for each of the configurations are summarized in the table below. 


Configuration Summary - Gear Up 


Config. 

Number 

Weight 

(lbs) 

Altitude 
(Ft MSL) 

Airspeed 

(Knots) 

Flaps 

(%) 

CG 

% MAC 

1 

140,000 

4,000 

1 60 

0 

20.85 

2 

140,000 

4,000 

145 

30 

20.85 

3 

160,000 

4,000 

1 75 

1 0 

20.85 

4 

1 40,000 

4,000 

1 55 

3 0 

f 20.85 


The transfer functions were obtained from the quadruples using System Technology's CC Pro- 
gram 7 . These aircraft transfer functions are listed here with each respective row of numbers designat- 
ing the corresponding configuration transfer function values. The nominal configuration, number 1, is 
represented by values in each row 1 below. 

_ kjq(<k»/Kic) , i 
Ztlh*) “ * ^, ong 

ZUb.) * ^Jong 


3.07E-04 (0) (-I.I7E-05) (0.40) (.61) Conf#l 

_ 3.’oYe' ’ 04 " "(0)’ “a 36E - 06) ’ '(b"797,"0. 568) ’ Conf #2 

/l,N> 27571 E -04 " TO)’ ’ ’ 167292" (a644) ToiTf # 3 

2~59E - 04 ' ’(b)" ' T 2ME- 06) * ~(0 8 lT.~0.508) ’ Conf # 4 


3.63E-05 

(0) 

(0.203) . 

(0.370, 3 008) 

Conf #1 

2.36E-05 

(0) 

(0.303) 

(0.476,3.51) 

Conf #2 

~ 2.77E-05 

(0.167) 

(0.351 , 3.038) 


Conf #3 

1.9 JE - 05 

(0) 

( 1) 

(0.460 , 3.426) 

Conf #4 


(1.438E - 05) 

(3.918E -02,0.130) 

(0.652,1.382) 

Conf #1 

(I.101E -05) 

(7.423E- 02,0. 147) 

(0.596,1.375) 

Conf #2 

(3.949E- 02,0. 118) 

(0.649,1.301) 


Conf #3 

(1.878E-05) 

(7.190E -02,0.1 38) 

(0.588,1.279) 

Conf #4 
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Application of QFT to the Problem of Failed In-Flight 
Controllers During Approach and Landing of a B-720 Aircraft 
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Cal Poly State University 
San Luis Obispo, CA 93407 


Abstract 


Previous studies by NASA Dryden have shown 
the use of throttles for emergency flight control to be 
extremely difficult, especially for landing. Flight control 
using only the throttles to achieve safe landing for a large 
jet transport airplane, the Boeing 720, is investigated using 
Quantitative Feedback Theory. Results are compared to an 
augmented control developed in a previous study. The 
controller corrected unsatisfactory open-loop 
characteristics by increasing system bandwidth and 
damping, but improving the control bandwidth 
substantially proved very difficult. The pitch controller is 
robust in conditions of no or moderate turbulence. The 
roll controller performed well in conditions of no 
turbulence, but is sensitive to moderate turbulence. 
Handling qualities of the augmented control for approach 
and landing were evaluated by piloted simulation flights. 

Notation 

TOFC Throttle-Only Flight Control 

QFT Quantitative Feedback Theory 

Cmu Velocity-pitch rate derivative 

Qp Sideslip-roll coupling derivative 

q pitch rate (deg/sec) 

y flight path angle (deg) 

P angle of sideslip (deg) 

cp bank angle (deg) 

6 pitch angle (deg) 

z thrust (lbs) 

CO i natural frequency 

^ damping ratio 

&Tc stick input(ful] defiection=l unit) 

pitch rate feedback gain 
flight path angle feedback gain 
sideslip angle feedback gain 
bank angle feedback gain 
G° n ut transfer functions 


(a) short form of (s+a) 

I£,a>] short form for s ^ + 2 £cos + co ^ 

c.g. center of gravity 

Introduction 

Through throttle manipulations, engine thrust has 
been found useful in providing some controllability for 
multiengine aircraft in emergency situations with severe or 
complete flight control system failures. This paper focuses 
on a particular application, a simulation augmented control 
developed by NASA for Boeing-720 aircraft Throttles-Only 
Right Control (TOFC) ^ NASA has found the use of 
throttles feasibl<rfor emergency flight control for a range of 
airplanes, and their analyses for a variety of aircraft are 
available in the literature ^“6. This controller was 
implemented on a high fidelity B720 flight simulator and 
obtained generally good pilot ratings by increasing the bare 
airframe Dutch-roll and phugoid damping. The primary 
aim of this study is to present an alternative control design 
technique based on Quantitative Feedback Theory(QFT) to 
further improve the Dutch-roll damping and to increase 
controller bandwidth for better handling qualities. 

The QFT technique ^ was chosen because of the 
insights it provides throughout its design process. It allows 
designers to specify a desired closed-loop frequency 
response with a desired control bandwidth and damping 
characteristic. It shows why the desired performance may 
not be achieved within the given control actuation and rate 
limits. 

In this paper a QFT augmented throttle-only flight 
path controller for approach and landing is presented. 
Complete details of the aircraft model and justification for 
TOFC are not included, but the reader is reminded that 
"trimming" must be possible and "controllability" must 
exist. Augmented control design using QFT is presented in a 
summary fashion. The full justification and step by step 
procedure may be found in Reference 8. 
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B-720 Linear Model 


Performance Sp eci fic ation 


The empirical transfer function developed for the 
engines is given in short form notation by 

G i(lbs) 275 

(.55X5) 

Severe bandwidth attenuation occurs beyond frequencies of 
1 rad/sec. For this application his prevented the increase the 
closed-loop bandwidth beyond 1 rad /sec within the range of 
available thrust (see Ref. 6). 

Four configuration variations for the B-720 were 
considered as described in the Appendix. They are 
characterized in both the longitudinal and lateral axes by 
excessive resonance, low phase and gain margins, low 
crossover frequency, and large phase angle roll-off. 

QFT Controller Design 

To apply QFT, the aircraft model is rearranged in a 
unit feedback form as shown in Figures 1 and 2. The inner 
pitch rate and sideslip loops were closed using Kq=60 and 
Kp = 4, which were the settings for the original simulation 
augmentation scheme. 


QFT allows designers to specify a desired closed- 
lo ^p frequency response with an upper bound Bu a lower 

bound B^and a tolerance <5 R specified to obtain robust 
performance. The maximum M m is also given to obtain a 
desired system damping. 


Table 1. QFT Performance Specification 


Freq.(r/s) 

0.1 

0.3 

05 

0.7 

1.0 

2.0 

5.0 

Bu(dB) 

17.0 

17.0 

17.3 

-16.0 

-4.0 

20 

-13.0 

BL(dB) 

16.8 

-15.0 

12.3 

4.6 

-12.4 

7.1 

-23.0 

S„ (dB) 

0.2 

2 

5 

8.4 

8.4 

9~1 

15.0 


The performance specification shown in Table 1 are 
the desired closed-loop responses for both the y - and <p - 
loops. These two feedback loops are piloted open-loop 
systems. Additional specifications are ususallly given for 
piloted systems, such as a desired control bandwidth of 2 
rad/sec. (see Ref. 9 for transport aircraft landing 
requirements) and a k/s slope near the crossover frequency. 
These added requirements promote good pilot handling 
qualities. 10 
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Design Constraints 


Four configurations were used to study the 
approach and landing of B-720 throttle-only flight control as 
summarized in Table 2 and in the Appendix. Configuration 
1 was used as the nominal confiuration for control design. 


Table 2, Flight Configurations for B-720 
(Gear Up) . 


Conf 

Weight 

(lbs) 

Alt 

(Ft) 

(MSL) 

A/S 

(Kts) 

Flaps 

(%) 

1 

140,000 

4,000 

160 

0 

2 

140,000 

4,000 

145 

30 

3 

160,000 

4,000 

‘ 175 

0 

4 

140,000 

4,000 | 

155 

30 


The performance bounds constraint is a curve on 
the Nichols Chart that shows the performance tolerance, 
from Table 1 at each specified frequency. Satisfying the 
tolerance constraint guarantees that the variation of the 
system response due to plant uncertainties will be no greater 
than <5 r There is a performance bound for each frequency. 

U contours are also shown on the Nichols Chart. 
The U contour is a M -circle that has the magnitude of M m , 
with part of the circle stretched for uncertainty at high 
frequencies. -By having the open-loop frequency response not 
penetrate the U contour, the system's damping will be 
guaranteed no less than the damping selected for M m as a 
design constraint. 


The of the nominal configuration(config. 1) is: 




(.562) [.624, .111] [.441, 1.57] (5.25) 


and the min. and max. Gg (d ( ^ are : 

■ _ .0053 (.162) [.35 ,3.011 

®- <dc£) ‘ (.40) [.42 , 1.48] [.66 , .01] (5.19) 

C r<dcg) -01 (-28) [.46, 3.43] 

e.(dcfi) • j 45 _ j_ 57 j [ 92 , 14 ] ( 5 24) 

Figure 3. Longitudinal Mode Parameter Variation 


The ^^cg) of the nominal configuration(config. 1) is: 


.09 [.47 , 3.65] 


C A.Wcg) notninal [ gl ( 15 j [ 26 _ i.07] (5.02) 

and the min. and max. of . are : 

P..(dcg) 

^♦(deg) .06 [.45 , 3.65] 

A.(d=g) mm _ (.98) (.60 , .15] [.24 , .93] (5.01) 

G * ( d eg) max -09 [-61 , 4.331 

/U<J»8) max ' (i.03) [1.0 , .20] [.29 , 1.09] (5.02) 


Figure 4 Lateral Mode Parameter Variation 


For inner-loop transfer functions G£ and G! the 

An An ' 

parameter variation given by the four configurations can be 
expressed minimums and maximums. There are tradeoffs 
between plant parameter variations and performance. The 
wider the parameter variation, the more restricted the 
constraints; consequently, more compensation is required. 

In this application, due to the engine response limits, the 
performance specification will have to be relaxed because 
there is not enough control power to provide all the 
compensation that is required to meet specifications. 

The minimum and maximum values of transfer 

functions, and G^ , form the uncertainty template. 

The QFT control package ^ allows the designer to input 
maximum and minimum plant parameter variations, but 
due to the software's limitation of handling the quantity of 
uncertainties, some of them were averaged. These variations 
are listed in Figures 3 and 4. 


Poles/zeros/gain compensation will be required to 
reshape the open- loop transfer functions from G£ and 

G* . On the Nichols Chart, gain raises the transfer function 

curve, a zero bends the curve to the right, and a pole bends 
the curve to the left. The compensation -selected forms the 

controller, G c . G^ and G^ after reshaping, become 

respectively and L^, the open-loop transfer functions 
of the flight path angle and bank angle feedback loops, 

wtae 4, = G£ * and I*, = Gj; * G'_ . L' , 

and should be kept on and above the Bo(jcaj), for each 
frequency, co; , on U er and L* 0 to assure robust 

performance. U eyr and must also not penetrate the U 

contour in order to obtain the desired damping. In this 
application the additional constraint existed which required 
the controller to be physically realizable (zeros not 
outnumbering poles). 
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Long i tudinal F li ght Pa t h Angle Controller 


Lateral Bank Angle Controller 


Transfer function G^ and its performance bounds, 
Bo(jcoj), and U contour are displayed on a Nichols Chart in 
FigureS. Since all the frequency points for G£ are below 

In 

their corresponding Bo(jo>j).reshaping is required. Pure gain 
compensation First raises the curve until it touches the U 

contour as shown in Figure 6(G^ ,n =16). Compensation is 
then added to avoid the U contour while satisfying all the 
Bo(j(Oj) constraints. Note that no realizable polcs/zero 
compensation could be found to do this. A zero at .1 rad /sec, 

for example, pulls the whole G^ curve to the right of the U 
contour. But then a pole at any location would make the 
G^ curve penetrate the U contour. Since a compensator 
with only one zero is physically ununrealizable for this 
application, the compensator, G c , for the flight path angle 

feedback loop, G e ‘" , is a pure gain of 16. The performance 
bound, hence the system robustness, was left unsatisfied. 

The frequency response of the close-loop transfer 
function, TJ , where 

Mn 

T i = L r ., / (1 + L' r ) = (G£ * G'J / (1 + g;- * G'_ ) 

is shown in Figure 7. It can be seen in Figure 7 that 5 Jf the 
spread between T max and T m j n , has exceeded the 5 0 over 

the frequency range .1 to .7 rad/sec. This is becauseLI did 

not satisfy the performance bounds over that frequency 
range. The frequency plot of the close-loop response after 
adding a prefilter is shown in Figure 8. A pure gain prefjltcr 
of 6 proved most effective in increasing the bandwidth and 
in meeting the prescribed specification. 


Transfer function G^ and its performance bounds 
Bo(ja>i), and U contour were similarly analyzed on a Nichols 


Chart. G* was not only below all the Bo(je)j) but also 
penetrated the U contour. A controller, 

Gf* =(s-k15)/(s+15), was added to G^ to reshape it. 
is shown on a Nichols Chart in Figure 9 and t..„ requency 
plot of the close-loop transfer function, T * , where 


T* = L f /(l + L* ) = fG P|n *G 9 ) is 

*i n eo ' v > , ' u e«i °/3 |n ) > is 

shown in Figure 10 which show's the close-loop frequency 


response of T* after reshaping but with no prefilter applied 

yet. Sufficient gain is available here. A lead compensator of 
(S+l)/(S-*-2) is added to haunch up the severely deteriorated 
curve at frequencies over 1 rad /sec, and a lag compensator 
of (S+.25)/(S+.15) is added to steeper the gain curve at low' 
frequencies to provide a smoother k/s curve for good pilot 
handling qualities. The close-loop response after adding the 
prefilter is similar to Figure 8 , The prefilter selected is 
15(S+25)(S+l)/((S*.l5)(S+2)>. 


Results and Discussion 


For bank angle control using only the throttles, p 
feedback was found effective in increasing Dutch-roll 
damping. Bank angle feedback is crucial to lateral phugoid 
damping. Yaw rate feedback, on the other hand, increases 
Dutch-roll damping very little and actually decreases the 
lateral phugoid damping. Tables 4 and 5 compare the 
dynamic modes of the bare airframe with those from the 
augmented control implemented in a previous simulation 
study (References 1 and 2 ), and with those from the QFT 
implementation. 


Table 3. Longitudinal Mode Comparison 



Density 

Phugoid 

Short 

Period 

Engine 

f 

Pre- 

Filter 


G Sj ‘ 

K r 

K , 

Bare 

Airframe 

(1.4E-6) 

(.04,. 13) 

(.65,1.4) 

(.55X5.2) 

* 

* 

* 

X 

X 

Simulation 

Augmented 

Control 

(4.7E-6) 

(.52,. 24) 

(.52,1.5) 

(.4)(5.2) 

10 

1 

10 

1 

4 

QFT 

Augmented 

Control 

(3.4E-6) 

(.62,32) 

(.46,1.6) 

(3X5.2) 

6 

16 

1 

1 

60 
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Table 4 Lateral Mode Comparison 



Spiral 

Dutch 

Roll 

Roll 

Engine 

Pre- 

Filter 

G i 

' G < 

K c 

Kp 

k p 

Bare 

Airframe 

(1.1E-4) 

(.12,-99) 

(1) 

(.55X5) 

* 

x 

X 

X 

X 

X 

Simulation 

Augmented 

Control 

(.73,35) 

(.15,-99) 

(1) 

(5) 

40 

l 

1 

.5 

1 

.5 

QFT 

Augmented 

Control 

(39) 

! 

(.29,1.0) 

(13) 

(.45X5) 

2.5C25) 

(1.25) 

l 

1 

(.15) 

(1.5) 

4 

X 


For longitudinal control, pure gain compensation 
was used. The short period mode has a frequency near 1 5 
rad /sec, which is beyond the frequency that the throttles 
can control. Therefore, the primary concern was to increase 
phugoid damping and frequency. The short period 
damping decreased from .52 to .46 rad/sec while the 
phugoid damping and frequency increased from .52 to .62 
and from .24 rad /sec to .32 rad /sec, respectively. This 
increase of response frequency can also be easily seen on the 
flight path angle response (not shown here). 

For lateral control where pole /zero compensation 
was used, the Dutch-roll damping was almost doubled, 
from .15 to .29. The "simulation augmented controller" 
caused a lateral phugoid mode, (.73, .35), which combines 
the spiral and the slow engine mode. The QFT controller 
eliminated the lateral phugoid mode and resulted in higher 
damping for the Dutch roll mode (0.29 versus 0.15). 

Turbulence Response 

The response of the QFT flight path controller 
under intermediate turbulence was excellent. Since gusts 
were input, more than one simulation run was made to 
examine the tracking integrity under turbulence. Both 
controllers performed well, but the lateral QFT controller 
showed undue sensitivity to Kp during bank angle tracking. 

During investigation of the bank angle tracking 
problem, it was found that the (3 being feedback in the 
simulation was the J3 at the c.g. instead of at the nose 
boom. The nose boom p was then modeled into the B-720 
simulator, which improved the bank angle tracking under 
turbulence. The P at the nose boom has two extra terms, one 
a function of roll rate, the other a function of yaw rate. It 
was thought that the extra yaw rate term might have 


stabilized the bank angle tracking (the roll rate term was too 
small and was neglected.) 

System response to configuration variations, for 
flight path control and for bank angle control, was very 
good. The robustness of the flight path control was 
improved by QFT. The Dutch-roll oscillation in the original 
simulation compensation was taken out by QFT 
compensation; however, the tracking did not improve. 

Conclusions and Recommendations 

For throttles-only pitch control using a QFT 
controller, the control bandwidth , tracking and control 
robustness were improved by QFT. For bank angle 
control, QFT has improved the Dutch-roll oscillation 
problem and performed well under no turbulence. 

However, the lateral phugoid tracking under intermediate 
turbulence did not perform well. Apparently a compromise 
is required between Dutch-roll and lateral phugoid mode 
damping. Further investigation is recommended for bank 
angle tracking under turbulence. The impact of system non- 
linearities, such as rate and thrust limits, was significant and 
resulted in a decrease in the bandwidth specification used in 
the QFT analysis. 
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Longitudinal Transfer Functions; 


r <?(deg/ sccj 
i( lbs) 


z(lbs) 


l ^ tong 


G 


y(dcg) 

z(lbs) 


Z(tb5) 


f &long 



236E-04 

(0) 

H.17E- 05) 

vu. .0) 

(0.61) 

1 sec) 233E - 04 

<D) 

C1-4H- 06) 

(0.635, 

0.563) 

z(!bs) 

1.976E- 04 
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